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SUMMARY 


The  objective  of  Task  1  of  this  program  is  to  evaluate  and 
characterize  laser- induced  processing  of  semiconductor  materials 
and  devices.  The  work  described  here  was  performed  under  Contract 
No.  MDA  903-78-C-0284  during  the  period  of  February  1979  to 
July  1979.  The  previous  work  performed  during  the  period  of 
June  1978  to  January  1979  was  already  described  in  the  semi¬ 
annual  technical  report  "Laser  Assisted  Semiconductor  Device 
Processing,"  No.  33292.000. 

The  TRW  contribution  to  laser-induced  processing  of  semi¬ 
conductor  surfaces  consisted  of  the  demonstration  of  the  effect 

of  lattice  damage  on  the  laser- induced  heating,1  and  the  develop- 

2 

ment  of  annealing  techniques  for  lightly  damaged  materials. 

This  new  technique,  pioneered  at  TRW,  consists  of  combining  a 
laser  anneal  with  a  short  low  temperature  thermal  anneal.  It 
exploits  the  unique  ability  of  laser  processing  to  dissolve 
lattice  defects  without  providing  sufficient  energy  for  melting 
or  regrowth.  A  detailed  discussion  of  the  effect  of  damage  on 
the  laser- induced  heating  is  given  in  Section  2.1  and  Appendix  A. 

In  our  new  technique,  the  sample  is  first  irradiated  by  a 
pulsed  laser.  This  removes  lattice  defects.  This  defect-free 
material  is  then  subjected  to  a  moderate  thermal  anneal.  This 
results  in  full  electrical  activation  without  significant  change 
in  the  impurity  distribution.  Thus,  by  combining  laser  irradiation 
with  short  low  temperature  thermal  anneal,  it  is  possible  to  both 
eliminate  defects  and  electrically  activate  lightly  damaged  semi¬ 
conductors.  Since  melting  is  avoided,  the  dopant  distribution 
remains  as  implanted.  The  results  of  these  experiments  are 
summarized  in  Section  2.2  and  Appendix  B.  Preliminary  experi¬ 
ments  have  also  shown  that  this  technique  can  be  applied  to  high 
fluence  implants. 

The  objective  of  Task  2  is  to  evaluate  the  laser  blow-off 
ion  implantation  source. 


* 


This  new  ion  implantation  technique  utilizes  pulsed  laser 
irradiation  of  a  solid  surface  of  pure  dopant  material  to  pro¬ 
duce  a  high  temperature  plasma  of  dopant  ions.  The  plasma  ex¬ 
pands  into  vacuum,  reducing  the  charge  density  to  the  point  where 
an  externally  applied  electric  field  can  be  used  to  extract  ions 
from  the  plasma.  In  principle,  these  ions  can  be  accelerated  to 
an  arbitrarily  high  energy  and  can  be  used  to  implant  semicon¬ 
ductor  materials.  A  subtle  aspect  of  this  method  is  the  pulsed 
nature  of  the  ion  beam  current.  For  reasonable  values  of  ion 
energy  and  ion  areal  density,  it  appears  that  annealing  temper¬ 
atures  can  be  reached  in  the  implanted  regions  for  very  short 
times  at  each  laser  pulse. 

Measurements  have  been  made  and  results  obtained  largely  in 
the  order  of  laser  performance,  implant  plasma  production,  plasma 
plume  characteristics,  ion  current  and  areal  uniformity,  implant 
material  deposition  without  acceleration,  and  implant  with  high 
energy  implant  ions. 
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TASK  1.  LASER  INDUCED  SURFACE  PROCESSING 


1.  INTRODUCTION 

In  our  work  at  TRW,  we  have  investigated  the  important 
difference  between  annealing  of  heavily  and  lightly  damaged 
materials.  Our  original  results  confirmed  findings  recently 
reported  in  the  literature:  for  lightly  damaged  samples,  the 
crystal  regrowth  was  poor,  and  only  partial  electrical  acti¬ 
vation  was  obtained,  though  the  displacement  damage  was  removed. 
In  contrast,  for  high  fluence  samples,  excellent  activation  was 
obtained . 

The  mechanism  for  the  Q-switched  laser  annealing  is  believed 
to  involve  melting  of  the  amorphized  surface  followed  by  liquid 
phase  regrowth.5’4  Here,  the  defects  are  almost  completely  re¬ 
moved  and  a  full  activation  of  dopants  at  concentrations  far 
exceeding  the  equilibrium  limit  of  solid  solubility  has  been 
reported.5  This  method  has  its  biggest  potential  in  applications 
where  slight  changes  in  the  implanted  profiles  can  be  tolerated. 
This  is  because  melting  results  in  some  impurity  redistribution. 
For  ruby  annealed  silicon  samples,  typical  profile  widths  are 
around  0.3  pm.*’5 

Our  early  experiments  lead  us  to  thoroughly  investigate  the 
effect  of  damage  on  the  laser- induced  heating.  For  the  theoreti¬ 
cal  effort,  we  have  developed  a  model*  which  explains  the  ob¬ 
served  differences  between  annealing  of  heavily  and  lightly 
damaged  surfaces.  This  is  done  by  taking  into  account  diffusion 
of  the  photoexcited  carriers.  This  is  discussed  in  Section  2.1 
and  Appendix  A.  In  our  experimental  effort,  we  have  subsequently 
developed  a  novel  technique  which  permits  laser  annealing  of 
lightly  damaged  materials.  This  has  the  very  attractive  feature 
of  maintaining  the  as-implanted  impurity  profiles.  This  is 
discussed  in  Section  2.2  and  Appendix  B.  Preliminary  experiments 
have  also  shown  that  this  technique  can  be  applied  to  high 
fluence  implants. 


2 .  RESULTS 


2.1  THE  EFFECT  OF  DAMAGE  ON  LASER  HEATING 

The  conversion  of  optical  to  thermal  energy  during  laser¬ 
annealing  consists  of  two  major  processes: 

•  Absorption  of  laser  light 

•  Conversion  of  absorbed  optical  energy  to  heat. 

In  strongly  absorbing  semiconductors,  the  dominant  absorption 
mechanism  at  frequencies  higher  than  the  bandgap  frequency  is 
interband  transitions.  The  absorbed  energy  is  released  to  heat 
through  the  thermalization  of  the  photoexcited  carriers,  i.e., 
heat  is  distributed  into  a  depth  comparable  to  carrier  diffusion 
length.  In  good  quality  materials,  the  carrier  diffusion  lengths 
are  longer  than  the  optical  absorption  depth.  Thus,  the  heating 
volume  is  enlarged,  and  for  the  same  laser  power,  the  average 
temperature  is  lower  than  in  heavily  damaged  material. 

The  carriers  recombine  preferentially  at  defect  sites, 
producing  localized  thermal  spikes.  If  the  density  of  the  non- 
radiative  recombination  centers  is  high,  the  average  temperature 
is  high.  This  is  the  case  in  heavily  damaged  materials.  At 
these  elevated  temperatures  epitaxial  regrowth  (liquid  phase) 
and  impurity  activation  occur.  If  the  density  is  low,  the 
average  temperature  is  low.  This  is  the  case  in  lightly  damaged 
materials.  Here  the  application  of  laser  fluxes  equivalent  to 
that  used  in  high  fluence  implants  produces  no  melting  in  low 
fluence  implants.  No  regrowth  and  little  or  no  dopant  acti¬ 
vation  occur  even  though  the  defects  dissolve  due  to  the  pref¬ 
erential  supply  of  energy  at  the  defect  sites.  Thus,  it  is 
possible  in  lightly  damaged  materials  to  remove  the  displacement 
damage  even  without  providing  sufficient  average  energy  to 
activate  implanted  impurities. 

Our  detailed  thermal  analysis  supports  this  picture.  We 
have  developed  a  model*  which  explains  the  observed  differences 
between  annealing  of  lightly  and  heavily  damaged  materials  by 


> 


i 

E 


a  w 

a 


i 


» 


taking  into  account  the  diffusion  of  the  photoexcited  carriers. 
Previous  thermal  models  were  concerned  with  heavily  damaged 
materials,  in  that  they  assumed  that  the  absorbed  optical  energy 
is  converted  to  heat  in  the  immediate  vicinity  of  the  absorbed 
photons.  In  our  work,  the  time  evolution  of  temperature  pro¬ 
files  was  calculated  for  pulsed  ruby  and  Nd:YAG  illumination  of 
silicon,  and  ruby  illumination  of  GaAs.  The  material,  optical, 
and  electronic  properties  were  varied,  to  represent  varied 
material  quality.  Surface  recombination  of  the  carriers  was 
also  taken  into  account.  The  effect  of  the  damage  on  the 
heating  was  found  to  be  much  more  pronounced  for  Nd:YAG  than 
for  ruby  irradiation  because  a  larger  portion  of  the  ruby  laser 
energy  (excess  above  bandgap)  is  given  rapidly  to  the  lattice 
through  quasi-thermalization  inside  the  semiconductor  energy 
bands.  Thermal  spikes  were  formed  as  the  result  of  the  prefer¬ 
ential  recombination  of  the  carriers  at  the  surface.  For  a 
sufficiently  fast  surface  recombination,  the  temperature  profile 
may  develop  a  secondary  maximum  in  the  bulk. 

2.2  LASER  ANNEALING  OF  LOW  FLUENCE  ION- IMPLANTED  SILICON 

The  most  successful  applications  of  Q-switched  lasers  have 
been  for  annealing  of  high  fluence  implants.^  The  method  involves 
a  pulsed  laser  such  as  a  Q-switched  Nd:YAG  or  ruby  laser,  with 
a  typical  pulse  length  between  10  and  200  nsec,  and  energy 
density  between  .5  and  3  J/cm  .  The  annealing  mechanism  is 
believed  to  be  melting  of  the  amorphized  surface  followed  by 
liquid  phase  regrowth. Here  there  is  almost  complete  removal 
of  defects  and  a  full  activation  of  dopant  concentrations  far 
exceeding  the  limit  of  solid  solubi 1 i ty . ^ ^  However,  this  method 
can  be  only  applied  where  slight  changes  in  implanted  profile 
can  be  tolerated.  This  is  because  melting  results  in  some  im¬ 
purity  redistribution.  Typical  profile  widths  in  ruby  laser 
annealed  silicon  are  around  0.3  ym.5,6 

The  second  laser  anneal  method  employs  a  cw  laser  operated 
in  a  scanning  mode.  For  example,  a  7W  cw  argon  laser  beam 
scanned  at  2.76  cm/sec  and  producing  lines  of  22  ym  width  was 
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used  to  anne&l  to  a  1000  X  depth  in  as- implanted  silicon.7 
The  recrystallization  is  through  solid  state  regrowth,  because 
of  the  low  powers  and  comparatively  long  anneal  times.  The 
dopant  profiles  are  not  changed  by  cw  laser  annealing. 

The  as- implanted  impurity  profiles  can  be  maintained  even 

7 

in  pulsed  laser  annealing  by  using  techniques  developed  at  TRW. 
These  techniques  consist  of  combined  laser  irradiation  and  low 
temperature  thermal  anneal.  In  our  work  at  TRW  we  have  system¬ 
atically  investigated  the  differences  between  annealing  of 
lightly  damaged  materials  and  heavily  damaged  materials.  We 

_  7 

have  applied  ruby  laser  pulses  of  1.5  to  2.5  J  cm  to  silicon 
surfaces  oriented  in  [100]  and  [ill]  directions  and  implanted 
with  1014  to  1015  B  cm  ^  and  10^4  to  5  x  10*4  B  cm  7.  Boron 
implantations  were  at  50  keV,  phosphorus  at  80  keV.  It  was 
found  that  laser  energy  densities  used  for  high  fluence  implants 
did  not  produce  melting  or  regrowth  in  the  low  fluence  implants. 
Although  the  electrical  activation  was  poor  the  lattice  defects 
(e.g.,  dislocation  loops)  were  removed.  We  have  concluded  that 
the  average  temperature  induced  in  the  lightly  damaged  material 
was  substantially  lower  than  the  temperature  induced  in  the 
heavily  damaged  material.  This  permitted  only  a  small  fraction 
of  the  implanted  ions  to  migrate  to  substitutional  positions. 

The  damage  was  removed  through  the  preferential  heating  at  the 
defect  sites. *  However,  when  a  sequential  moderate  thermal 
anneal,  10  min  at  1000°C,  was  applied  to  these  samples,  a  full 
electrical  activation,  and  without  significant  impurity  redis¬ 
tribution,  was  obtained.  Thus,  by  combining  a  laser  anneal 
with  a  short  thermal  anneal,  it  was  possible  to  achieve  defect 

elimination  as  well  as  complete  electrical  activation.  This 

? 

was  the  first  reported  successful  use  of  pulsed  laser  radiation 
for  annealing  of  low  fluence  implants. 

Preliminary  experiments  have  shown  that  we  can  also  apply 
these  techniques  to  high  fluence  implants.  One  type  of  im¬ 
plantation,  that  of  5  x  1014  P  cm  7  into  a  [111]  surface, 
resulted  in  amorphism.  Our  early  experiments  as  well  as  others 
in  the  literature  showed  that  irradiation  of  amorphized  surfaces 


by  a  Q-switched  ruby  laser  results  in  melting.  Thus,  in  order 
to  avoid  melting,  the  damaged  region  should  first  be  recrystal¬ 
lized.  This  was  done  by  introducing  a  moderate  thermal  anneal 
before  the  irradiation.  To  date,  the  best  results  were  obtained 
with  a  two-step  thermal  anneal,  30  min  at  550°C,  and  10  min  at 
1000°C.  This  caused  the  amorphous  layer  to  regrow  into  single 
crystallinity  by  solid  phase  epitaxy.  The  subsequent  irradi¬ 
ation  eliminated  defects  without  melting.  Not  only  did  we 
eliminate  the  defects,  but  by  avoiding  melting,  we  also  main¬ 
tained  the  desirable  as-implanted  impurity  profiles. 

The  demonstrated  techniques  have  an  enormous  potential  for 
practical  applications.  Low  fluence  implants  account  for  more 
than  50  percent  of  the  potential  implant  uses  since  all  bases 
for  junction  transistors  fall  into  this  category.  The  ability 
to  anneal  high  fluence  implants  without  dopant  redistribution 
is  important  for  the  production  of  shallow  emitters,  which  are 
essential  for  an  entire  category  of  devices,  such  as  microwave 
transistors.  The  ability  to  anneal  material  without  changing 
the  implanted  profiles  would  permit  full  exploitation  of  the 
ion-implant  unique  capability  of  creating  narrow  accurately 
controlled  impurity  profiles. 


TASK  2.  LASER  BLOW-OFF  ION  IMPLANTATION 


3.  INTRODUCTION 

A  new,  potentially  less  expensive,  and  relatively  simple 
technique  for  ion  implantation  of  extended  areas  of  semiconductor 
material  has  been  conceived  and  is  being  experimentally  developed 
under  contract  to  DARPA.  This  new  method  involves  a  source  of 
implant  ions  produced  by  laser  "blow-off"  from  a  piece  of  dopant 
material.  Energy  from  a  short-pulse  laser  is  focused  onto  the 
surface  of  the  dopant  material,  forming  a  small  bubble  of  dopant 
plasma  which  "blows  off."  This  plasma  bubble  is  allowed  to 
expand  into  vacuum  until  its  density  and  area  become  suitable 
for  electrostatic  acceleration  to  the  necessary  ion  implant 
energy.  Key  features  include  the  easy  formation  of  intense 
ion  beams  of  virtually  any  dopant  material,  with  beam  areas  of 
several  inch  diameter,  and  without  the  necessity  for  slower  and 
more  expensive  rastering  of  ion  beams  from  sophisticated  ion 
beam  machines.  Perhaps  most  significant,  although  also  more 
speculative,  is  the  potential  for  self-anneali ng  of  implant 
damage  during  the  implant  process.  This  possibility  arises 
from  the  unique  situation  of  the  short  pulse  of  the  high  inten¬ 
sity,  high  energy  ion  beam,  which  rapidly  deposits  energy  exactly 
in  the  implant  region  to  momentarily  raise  the  material  temper¬ 
ature.  If  this  possibility  can  be  realized  in  practice,  an 
entire  and  sometimes  troublesome  separate  annealing  process 
normally  associated  with  ion  implantation  can  be  eliminated. 

Thus  far  an  apparatus  has  been  built  and  tested  for  implanting 
ions  by  this  method.  Implants  of  boron  ions  into  silicon  substrates 
have  been  made  at  energies  up  to  50  keV.  The  upper  energy  was 
limited  by  electrical  breakdown  across  inadequate  vacuum  feed¬ 
throughs,  a  problem  which  is  being  corrected.  The  quality  of  the 
implants  is  being  assayed  with  Ion  Microprobe  Mass  Analysis  (IMMA) 
equipment.  The  IMMA  measurements  are  ideal  for  determining 


precisely  the  implant  ion  depth  profiles  and  surface  areal 
uniformity,  and  provide  meaningful  comparisons  between  the  new 
and  the  conventional  methods.  In  fact,  IMMA  has  been  more 
appropriate  as  a  diagnostic  than  actual  semiconductor  device 
fabrication  tests  to  this  time,  since  the  success  or  failure 
of  device  fabrication  may  say  little  about  large  and  small 
scale  uniformity  and  overall  depth  profiles  which  have  been 
our  initial  concern.  Major  emphasis  has  therefore  been  placed 
on  the  IMMA  evaluation  technique  up  to  now. 

In  the  next  contract  period,  major  emphasis  will  be  shifted 
from  IMMA  as  diagnostic  to  the  quality  of  actual  semiconductor 
devices  fabricated  with  the  laser  blow-off  implantation  tech¬ 
nique  and  to  C-V  type  measurements.  This  shift  is  particularly 
important  in  studying  the  self- annealing  feature  of  the  new 
method  since  IMMA  can  say  little  about  semiconductor  structure 
and  activation  properties.  And  in  the  final  analysis  the 
quality  of  the  fabricated  semiconductor  device  operation  is 
clearly  the  key  measure  of  success. 


4.  TECHNICAL  DISCUSSION 


4.1  TECHNIQUE 

Ion  implantation  of  dopants  in  semiconductor  materials  has 
become  a  standard  practice  in  the  manufacture  of  solid  state 
electronic  devices.  Ion  implantation  is  a  non-equilibrium  pro¬ 
cess  that  permits  a  degree  of  control  over  the  spatial  distribution 
of  dopant  atoms  in  the  host  crystal  that  cannot  be  achieved 
through  thermal  diffusion.  Dopant  species  with  low  diffusion 
rates  even  at  elevated  temperatures  can  often  be  handled  with 
relative  ease  through  ion  implantation.  However,  a  deleterious 
side  effect  of  ion  implantation  is  damage  to  the  host  crystal 
caused  by  the  high  energy  ions  necessary  to  penetrate  to  the 
desired  depths.  High  temperature  annealing  is  commonly  used  to 
repair  implant- induced  damage.  In  some  cases  the  annealing 
process  itself  causes  a  degradation  of  the  implanted  material. 

For  example,  high  temperature  annealing  of  implanted  GaAs  results 
in  migration  of  As  atoms  from  the  crystal.  In  large  scale  inte¬ 
grated  circuits,  heat  treating  may  alter  the  properties  of  the 
intrinsic  material  at  locations  remote  from  ion  implanted  sites. 

The  following  part  of  this  report  describes  progress  towards 
developing  a  new  ion  implantation  technique  which  shows  promise 
of  providing  a  sel f- anneal ing  feature  as  well  as  being  relatively 
simple,  fast,  with  large  areal  coverage  and  amenable  to  a  very 
wide  range  of  dopant  species. 

This  new  ion  implantation  technique  utilizes  pulsed  laser 
irradiation  of  a  solid  surface  of  pure  dopant  material  to  produce 
a  high  temperature  plasma  of  dopant  ions.  The  plasma  expands 
into  vacuum,  reducing  the  charge  density  to  the  point  where  an 
externally  applied  electric  field  can  be  used  to  extract  ions 
from  the  plasma.  In  principle  these  ions  can  be  accelerated  to 
an  arbitrarily  high  energy  and  can  be  used  to  implant  semiconductor 
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materials.  A  subtle  aspect  of  this  method  is  the  pulsed  nature 
of  the  ion  beam  current.  For  reasonable  values  of  ion  energy 
and  ion  areal  density,  it  appears  that  annealing  temperatures 
can  be  reached  in  the  implanted  regions  for  very  short  times  at 
each  laser  pulse.  Thus  the  possibility  presents  itself  that  the 
use  of  high  intensity  short  duration  ion  beam  pulses  for  ion 
implantation  may  be  accompanied  by  self -annealing  of  the  implant 
damage.  One  of  the  main  objectives  of  this  program  is  to  care¬ 
fully  examine  this  possibility  experimentally  by  systematically 
evaluating  the  characteristics  of  semiconductor  devices  implanted 
with  this  laser  blow-off  technique  over  a  variety  of  practical 
conditions . 


4.2  EXPERIMENTAL  APPARATUS 


A  schematic  diagram  and  picture  of  the  ion  implantation 
system  developed  here  is  shown  as  Figures  la  and  lb.  This 
system  is  an  adaptation  of  the  high  intensity  metal  ion  beam 
source  described  in  Reference  8.  The  ions  to  be  implanted  are 
created  by  pulsed  laser  irradiation  of  a  solid  block  of  the  im¬ 
plant  material.  In  the  present  case  a  Nd:YAG  laser  is  Q-switched 
with  a  Pockels  cell,  producing  a  1  joule  laser  pulse  with  less 
than  0.1  microsecond  length  at  a  wavelength  of  1.06  micron.  The 
laser  beam  passes  through  steering  prisms  and  vacuum  system  win¬ 
dows  and  is  focused  to  a  2  mm  diameter  spot  on  the  implant  material 
by  a  60  cm  focal  length  lens.  Tests  thus  far  have  utilized  a 
solid  boron  target.  This  laser  target  is  located  in  an  electric 
field-free  region,  with  the  normal  to  its  surface  directed  approxi¬ 
mately  toward  the  silicon  wafer  to  be  implanted.  The  laser  beam 
strikes  the  surface  from  an  angle  of  about  25°  to  the  normal. 

The  position  and  angle  of  the  boron  target  are  adjustable  to 
optimize  the  intensity  and  uniformity  of  the  ion  beam  in  the 
direction  of  the  silicon  wafer,  since  these  parameters  are  compli- 

O  Q 

cated  functions  of  the  laser  beam  energy  and  angle.  ' 

The  result  of  the  rapid  deposition  of  laser  energy  on  the 
implant  material  surface  is  to  produce  a  high  temperature  plasma 
"bubble"  of  implant  material  ions  and  electrons,  which  "blow  off" 
from  the  surface.  In  the  present  configuration  this  plasma  bubble 
expands  freely  until  it  reaches  the  ion  acceleration  region. 

This  distance  is  presently  37  cm.  In  the  first  measurements  no 
acceleration  was  utilized,  and  the  silicon  wafers  and  an  array 
of  current  measuring  Faraday  cups  were  placed  at  this  position. 

This  set-up  allowed  the  measurement  of  implant  ion  fluence  and 
areal  uniformity.  This  configuration  may  prove  to  be  valuable 
in  its  own  right  as  a  method  for  depositing  a  thin  layer  of  cer¬ 
tain  materials  for  subsequent  laser  diffusion  processing,  since 
atoms  of  essentially  any  solid  dopant  can  be  deposited  under 
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Figure  la.  Schematic  Diagram  of  Laser  Blow-Off  Ion  Implantation  System 


vacuum  conditions.  The  ion  energies  are  a  few  hundred  electron 
volts  in  this  case,  and  thus  barely  penetrate  the  surface. 

In  the  case  where  high  energy  implantation  is  to  be  used, 
acceleration  of  the  implant  ions  is  necessary  as  they  reach  the 
acceleration  region.  In  the  present  configuration,  this  accel¬ 
eration  occurs  between  the  grids  shown.  The  first  grid  is  held 
at  the  same  potential  as  the  field-free  plasma  expansion  region. 
The  second  is  held  at  a  high  negative  voltage  equal  to  the  ion 
kinetic  energy  in  eV  desired.  The  silicon  wafer  is  held  at  the 
same  potential  as  the  second  grid  at  a  distance  approximately 
one  centimeter  beyond.  This  "drift"  distance  is  intended  to 
wash  out  any  non-uniformity  in  implant  caused  by  grid  shadows 
or  focusing  effects.  The  drift  distance  cannot  be  too  long  be¬ 
cause  space  charge  spreading  becomes  severe  because  the  plasma 
electrons  are  rejected  by  the  ion  acceleration  region.  The  grids 
are  presently  100  mesh  stainless  steel. 

For  initial  tests,  two  types  of  silicon  wafers  were 
obtained  from  Monsanto.  Both  types  are  3  inches  in  diameter 
and  0.015  inches  thick,  are  oriented  100,  and  have  resistivities 
of  10  ohm- cent imeters .  One  type  is  N  doped  with  phosphorus 
and  the  other  type  is  P  doped  with  boron.  In  order  to  economize 
on  wafer  material,  the  wafers  are  cut  into  cm-size  chips  for  use 
when  the  measurements  permit. 

The  accelerating  voltage  is  supplied  by  a  Plastic  Capacitors 
Inc.,  50  kV  power  supply,  which  charges  an  0.05  microfarad  capac¬ 
itor  across  the  accelerating  grids.  In  the  event  of  arc  break¬ 
down,  the  grids  would  have  to  dissipate  60  joules  of  energy. 

2 

If  an  accelerating  region  area  of  100  cm  and  an  implant  ion 
density  of  1  microcoulomb/cm  /pulse  is  considered,  the  acceler¬ 
ating  voltage  drops  from  50  kV  to  48  kV  during  the  pulse.  A 
larger  capacitance  value  would  reduce  the  voltage  change,  but 
the  present  value  is  safer  during  development. 

The  vacuum  system  consists  of  fairly  standard  components 


including  a  10  inch  oil  diffusion  pump  with  water-cooled  baffle 
and  mechanical  forepump.  This  system  should  be  adequate  for 
developmental  purposes,  but  final  definitive  quantitative  device 
tests  may  necessitate  an  ultra-high  vacuum  system  to  minimize 
surface  contamination. 
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4.3  RESULTS 


Measurements  have  been  made  and  results  obtained 
largely  in  the  order  of  laser  performance,  implant  plasma  pro¬ 
duction,  plasma  plume  characteristics,  ion  current  and  areal 
uniformity,  implant  material  deposition  without  acceleration, 
and  implant  with  high  energy  implant  ions. 

The  Nd:YAG  laser  is  operated  at  an  output  energy  up  to 
1  joule  per  Q-switched  pulse.  This  is  lower  than  anticipated 
and  appears  to  involve  a  problem  with  the  polarizer  and  Pockels 
cell.  The  output  energy  also  varies  more  from  shot  to  shot  than 
had  been  expected,  about  -  5  to  10  percent.  This  leads  to  a 
roughly  similar  variation  in  implant  ion  fluence  from  shot  to 
shot.  Although  the  laser  energy  and  stability  are  less  than 
might  be  desired,  they  should  be  sufficient  to  carry  out  the 
present  tasks.  So  long  as  each  ion  pulse  is  monitored,  and  a 
number  of  shots  are  used  per  implant  dose,  it  should  be  possible 
to  control  the  total  implant  dose  within  the  necessary  limits. 

It  was  therefore  decided  to  accept  the  laser  performance  as  is 
for  now,  and  move  ahead  with  other  aspects.  It  is  certain  that 
better  laser  performance  is  possible,  and  that  eventually  when 
high  pulse  repetition  rates  are  needed  a  different  laser  will 
be  necessary  in  any  event. 

A  piece  of  pure,  solid  boron  of  dimension  of  order  1  cm 
was  obtained  for  the  laser  blow-off  target,  that  is,  the  source 
of  implant  ions.  The  actual  purity  is  thought  to  be  high,  but 
this  will  be  checked  mass  spectrometrical ly  during  analysis  of 
the  blow-off  deposited  material,  as  described  below.  Even  with 
many  laser  shots,  very  little  cratering  of  the  boron  occurred. 
This  is  consistent  with  earlier  work,  which  indicated  that  less 
than  one  micron  of  material  is  removed  per  shot.  No  apparent 
difficulty  was  encountered  in  the  implant  ion  plasma  production. 

In  order  to  examine  the  characteristics  of  the  expanding 


plasma  plume  with  regard  to  ion  energy,  current  and  number  den¬ 
sity,  and  angular  distribution,  an  array  of  5  Faraday  cups  was 
placed  symmetrically  along  an  arc  perpendicular  to  the  normal 
to  the  boron  laser  target  in  a  plane  containing  the  normal.  The 
cups  were  41  cm  from  the  target  and  were  spaced  4.2  apart  with 
respect  to  the  target,  that  is,  the  5  cups  observed  the  target 
at  angles  of  8.4°,  4.2°  and  0°  each  side  of  the  normal.  The 
distance  between  centers  of  the  outer  cups  was  12  cm.  The  Fara¬ 
day  cups  detected  ions  within  the  expanding  plasma  plume  as  the 
ions  reached  the  cup  collectors.  These  cups  were  designed  and 
tested  for  high  current  operation  at  the  plasma  densities  en- 
countered.  The  effective  detection  area  was  1  cm  .  These  cups 
were  of  planar  construction  so  that  the  ion  transit  time  within 
the  cup  was  negligibly  small  compared  to  the  transit  time  from 
the  laser  target  to  the  cup.  Elements  within  the  cup  were  electri¬ 
cally  biased  to  suppress  secondary  electron  emission.  The  Faraday 
cups  could  be  operated  in  either  the  current  sensing  mode  or  the 
total  charge  collection  mode,  with  the  signals  displayed  directly 
on  an  oscilloscope  in  either  case.  The  present  results  with 
boron  ions  are  consistent  with  earlier  theoretical  and  experi- 
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mental  work  done  in  this  laboratory  on  laser  blow-off  plumes. 

This  work  involved  a  model  based  on  the  following  physical  as¬ 
sumptions. 

It  was  assumed  that  the  deposited  laser  energy  vaporizes 
and  ionizes  Nq  particles  from  a  small  region  of  the  source  mate¬ 
rial  surface.  This  hot  plasma  "bubble"  of  Nq  ions  begins  to 
expand  against  the  substrate  and  imparts  momentum  to  it,  thereby 
giving  the  center  of  mass  of  the  N  ions  a  velocity  v  in  the 
laboratory  coordinate  frame.  Because  of  symmetry,  vcm  is  normal 

to  the  surface.  After  the  center  of  mass  of  the  N  ions  has 

o 

moved  a  distance  from  the  surface,  the  number  density  of  ions 
remains  high  enough  so  that  most  of  the  Nq  ions  equilibrate  at 
a  temperature  T  relative  to  their  center  of  mass,  which  continues 


I 

I 


i 
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to  move  at  v  .  That  is,  an  observer  moving  with  the  center  of 
cm 

mass  would  observe  an  essentially  spherically  symmetric  plasma 
"bubble"  in  equilibrium  at  T,  and  freely  expanding.  As  the  number 
density  continues  to  decrease,  the  ion  velocity  distribution 
corresponding  to  equilibrium  at  temperature  T  becomes  "frozen" 
in  the  center-of -mass  system.  A  laboratory  observer  at  a  larger 
distance  from  the  source  laser  spot  observes  ions  arriving  with 
velocity  essentially  along  a  line  from  source  spot  to  observer, 
but  with  a  plume-like  angular  distribution  due  to  the  center-of- 
mass  velocity  vcm  superimposed  on  the  purely  thermal  velocity 
distribution. 


The  theory  thus  consists  simply  of  calculating  the  number 
density  p  and  current  density  j"  of  ions  to  be  observed  in  the 
laboratory  frame  resulting  from  Nq  ions  expanding  with  spherical 
symmetry  and  having  equilibrium  velocity  distribution  corresponding 
to  T  in  their  center-of -mass  system,  and  with  a  superimposed 
center-of -mass  velocity  v 

cm 


Figure  2  shows  the  coordinate  system  for  the  theoretical 
model.  The  present  adaptation  differs  in  that  a  solid  target 
is  struck  by  the  laser  beam  on  the  front  surface.  The  current 
density  of  ions  J  in  ions  per  unit  area  per  unit  time  is  given 

y  +r  oil  l  5 1 2 \ / 2  \3/2/  kT  \1/2  No  /-3A2\  /  . 

j(r,0,A,n)  =  \~TT-)(tt/  (~1^/  ~3  exP  (  — 2“  )  *  ( <t,n) 


x  exp 


H 


-32  A  -  2 
( 4>n ) 


where 


■  1/ 


(A  cos  0) 2  ♦  y- 


1-3(A  cos  0)  (<t>n)  /  4 

(A  cos  0) 


H 


‘max  ■  3/8  0  V'”/3kT 


and 


n 

A  VRMS 


t/t 


max 


vcm  =  X  ^SkT/m 


19 


Figure  2.  Geometric  Arrangement  for  the  Theoretical  Model  and  Experiment 
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gives  the  relation  defining  the  "blow-off"  parameter  X.  Small  X 

implies  small  vcm  and  hence  a  broad  angular  distribution,  while 

large  X  implies  a  narrow  "plume".  The  other  parameters  above 

include  time  t,  ion  mass  m,  Boltzmann's  constant  k,  temperature 

T,  and  time  of  J  ov,  t  Figure  3  shows  these  results  graphi- 

cally  in  terms  of  J/J  versus  n,  and  therefore  represents  an 
/  +max  ’ 

expected  picture  of  j  versus  t  on  an  oscilloscope  trace  of  a 
current  sensing  Faraday  cup.  Picture  A  of  Figure  4  shows  an 
actual  trace  of  a  boron  ion  pulse  from  the  present  measurement. 

The  points  on  Figure  3  are  the  same  data  plotted  there.  The 

q 

skewness  away  from  the  theoretical  curves  has  been  noted  before 
for  ions,  while  the  agreement  is  quite  good  for  neutral  atom 
blow-off.^  The  X  values  depend  on  laser  energy  density  and 
the  target  materials,  and  are  thus  subject  to  some  control.  The 
general  agreement  with  theory  is  sufficiently  good  to  predict 
currents  and  densities  as  functions  of  time,  distance  and  angles. 

Of  particular  importance  is  the  sensitivity  of  the  angular  dis¬ 
tribution  to  X,  as  seen  next.  The  total  number  of  ions  passing 

/\ 

across  unit  area  normal  to  r  per  pulse  is  given  by 

/  Tdt  -  -2—2  /|  exp (-3X2/Z) 

4-nr 

x[X  +  /i(l/2  +  x2)  Cl  +  erf  x)  r 

X  =  /3/ 2  (X  cos  0) 

2_  X 

erf  x  =  f  exp(-u  )  du 

o 

Figure  5  shows  these  theoretical  angular  distributions  as  functions 
of  X.  Also  included  are  previous  experimental  points  obtained 
for  aluminum  ions  and  present  points  obtained  for  boron  ions. 


where 

and 
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values  for  boron  ions.  Boron  ion  current  vs  time 
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Figure  4.  Unaccelerated  Boron  Ion  Pulses 


Theoretical  Curves  for  J  3  dt  versus  Detector  Angle  0  at  Constant  Distance 
r  and  Normalized  to  0  the  Value  for  0=0.  Points  are  measured 
values  for  aluminum  and  boron  ions. 


It  is  obvious  that  small  values  of  A  are  necessary  if  one  is  to 
implant  large  areas  with  high  uniformity  without  "rastering" 
the  wafer  between  shots.  The  ion  arrival  time  profiles  may  also 
be  very  important  in  self -annealing  aspects,  since  the  thermal 
pulse  will  reflect  the  arrival  time  pulse.  Figures  4  and  5  from 
Reference  8  show  examples  of  widely  differing  angular  distributions 
depending  on  laser  beam  conditions.  The  presently  obtained  angu¬ 
lar  distributions  are  adequate  for  present  purposes  and  will  be 
optimized  further,  but  it  is  clear  that  variations  in  ion  fluence 
as  a  function  of  angle  put  eventual  limits  on  the  size  of  wafers 
which  can  be  implanted  with  high  uniformity  over  the  surface. 

Measurements  have  been  made  to  determine  the  actual  ion  cur¬ 
rent  and  areal  implant  ion  density  by  combining  Faraday  cup  array 
current  measurements  with  implant  ion  deposition  without  accel¬ 
eration.  The  currents  and  charge  collected  are  displayed  directly 
on  an  oscilloscope  for  the  two  outer  and  the  middle  Faraday  cups 
in  the  array.  The  apertures  of  the  other  two  cups  are  covered 
with  silicon  wafer  chips.  A  number  of  laser  shots  are  then  made 
to  deposit  boron  on  the  silicon  chips,  at  the  same  time  summing 
the  charge  collected  on  either  side  of  the  chips.  Table  1  shows 
the  results  of  one  such  run.  In  this  particular  case,  cup  A  was 
on  the  side  of  the  normal  toward  the  laser  beam,  and  the  results 
show  the  asymmetry  toward  the  laser  beam  which  was  reported 
earlier.  The  uniformity  of  +  2%  over  the  3  cm  interval  between 
cups  A  and  M  is  acceptable  for  now,  but  the  -  81  variation  between 
cups  B  and  M  is  not.  This  run  shows  a  particular  strong  depedence 
on  angle,  a  dependence  itself  dependent  on  laser  conditions  as 
noted  above. 

The  silicon  chips  deposited  with  boron  in  runs  like  Table  1 
were  analyzed  on  an  Applied  Research  Laboratories  IMMA  (Ion 
Microprobe  Mass  Analyzer).  A  schematic  diagram  of  this  instrument 
is  shown  as  Figure  6.  It  has  been  used  extensively  in  early  work 
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Figure  6.  Applied  Research  Laboratories 
(Ion  Microprobe  Mass  Analyzer) 


in  order  to  determine  implant  uniformity,  purity,  and  depth 
profiles.  Actual  semiconductor  device  implants  will  of  course 
be  necessary  to  prove  the  effectiveness  of  the  blow-off  implant 
method  in  the  final  analysis. 

Figure  7  shows  the  depth  profile  of  a  boron  implant  into  a 

3-inch  diameter  silicon  wafer,  with  an  implant  ion  energy  of 

40  KeV.  One  hundred  laser  shots  were  employed,  and  the  ion 

12  2 

fluence  was  a  few  times  10  boron  ions  per  cm  per  laser  shot. 
The  depth  profile  measurements  can  be  made  absolute  by  cali¬ 
brating  the  IMMA  against  a  known  profile  produced  by  conven¬ 
tional  ion  implantation  techniques. 
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APPENDIX  A 


THE  EFFECT  OF  CARRIER  DIFFUSION  ON 
LASER-HEATING  OF  LIGHTLY  DAMAGED  SEMICONDUCTORS* 

J.  Z.  WILCOX 

TRW  Defense  and  Space  Systems  Group,  Redondo  Beach,  California  90278 

ABSTRACT 

It  has  been  observed  that  the  surface  temperature  obtained 
during  laser  heating  of  lightly  damaged  material  is  generally 
considerably  lower  than  the  temperature  obtained  in  heavily  damaged 
material.  We  develop  a  model  which  explains  this  observation 
by  taking  into  account  the  diffusion  of  the  photoexcited 
carriers.  The  time  evolution  of  temperature  profiles  was 
calculated  for  pulsed  ruby  and  Nd:YAG  illumination  of  silicon, 
and  ruby  illumination  of  GaAs.  The  material,  optical,  and 

*  Supported  in  part  by  DARPA  Contract  MDA903-78-C-0284. 


electronic  properties  were  varied,  according  to  the  material  quality. 
Surface  recombination  of  the  carriers  was  also  taken  into  account.  For 
the  same  laser  power,  the  induced  temperature  was  lower  and  the  heating 
times  were  longer  in  good  quality  materials  than  in  heavily  damaged 
materials.  This  is  because  of  longer  carrier  diffusion  lengths  in  good 
quality  material.  The  effect  of  the  damage  on  the  heating  was  found  to 
be  much  more  pronounced  for  Nd:YAG  than  for  ruby  irradiation  because  a 
larger  portion  of  the  ruby  laser  energy  (excess  above  bandgap)  is  given 
rapidly  to  the  lattice  through  quasi-thermal ization  inside  the  semi¬ 
conductor  energy  bands.  Thermal  spikes  are  formed  as  the  result  of 
the  preferential  recombination  of  the  carriers  at  the  surface.  For 
a  sufficiently  fast  surface  recombination,  the  temperature  profile 
may  develop  a  secondary  maximum  in  the  material  interior. 


I.  INTRODUCTION 


There  is  considerable  interest  in  the  use  of  lasers  for  processing 
of  semiconductor  surfaces.^  The  majority  of  experiments  has  been  in 
pulsed  laser  annealing  of  comparatively  high-dose  implantations.  In 
these  experiments,  the  resulting  activated  dopant  concentrations  were 
often  in  excess  of  the  limits  of  equilibrium  solid  solubility,  and  the 

shapes  of  the  ion-implanted  profiles  were  changed  as  a  result  of  the 

2 

illumination.  This  led  to  the  suggestion  that  the  near-surface  regions 
of  the  irradiated  material  actually  melted  and  underwent  a  type  of  liquid 
epitaxy  in  these  experiments.  The  suggested  model  is  strongly  supported 
by  numerical  results  of  several  model  calculations. 

Heavily  damaged  material  has  been  successfully  annealed  also  with  scanned 
cw  Nd:YAG  lasers.^  No  significant  dopant  redistribution  was  observed 
after  the  annealing,  indicating  that  the  regrowth  was  through  solid  state 
epitaxy.  This  was  made  possible  because  of  the  long  times  involved  in 
these  experiments.  Nd:YAG  laser  annealing  was  found  to  be  a  very  sensi¬ 
tive  function  of  the  dose  of  the  implanted  impurities. 

Laser- induced  annealing  of  lightly  damaged  material  has  received  much 

c 

less  attention.  We  have  studied  the  use  of  Q-switched  ruby  lasers  for 

activation  of  low  fluence  non-amorphous  silicon  and  neutron-transmutated 

silicon,  and  found  that  powers  required  for  activation  of  non-amorphous 

surfaces  were  much  larger  than  those  required  for  amorphous  surfaces.  The 

required  powers  depended  strongly  on  the  dose  of  implanted  impurities. 
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For  the  neutron-transmutated  silicon,  an  energy  density  of  4  J/cm  was  not 
effective  for  activation,  unless  additional  damage  such  as  Ne  implantation 
was  introduced  into  the  material .  Lightly  implanted  wafers  showed  remelted 
spots  at  the  surface,  while  none  was  observed  at  surfaces  of  amorphized 


layers.  The  results  of  these  experiments  could  not  be  explained  in  terms 
of  a  weaker  optical  absorption  in  crystalline  silicon  alone. 

The  present  work  addresses  the  difference  between  laser-induced 
heating  of  heavily  and  lightly  damaged  materials.  It  presents  the  results 
of  the  thermal  calculations  which  take  into  account  the  diffusion  of  the 
photoexcited  carriers.  Previous  thermal  models  have  been  concerned 
with  heavily  damaged  materials,  in  that  they  assumed  that  the  absorbed 
optical  energy  is  converted  to  heat  in  the  immediate  vicinity  of  the 
absorbed  photons.  However,  in  strongly  absorbing  semiconductors,  the 
primary  effect  of  the  absorbed  radiation  is  to  generate  carriers.  Heat 
is  generated  when  the  carriers  recombine,  i.e.,  after  they  have  diffused 
a  distance  on  the  order  of  carrier  diffusion  length.  In  good  quality 
materials,  the  diffusion  lengths  may  be  longer  than  the  optical  absorption 
depth.  Hence,  the  carrier  diffusion  enlarges  the  heating  volume,  and,  for 
the  same  laser  power,  causes  the  surface  temperature  in  lightly  damaged 
material  to  be  lower  than  in  heavily  damaged  material.  Furthermore,  a 
large  number  of  carriers  will  reach,  and  recombine  at,  the  surface,  making  it 
qualitatively  feasible  to  melt  surface  defects  in  lightly  damaged  materials. 

The  conversion  of  optical  energy  to  heat  consists  of  two  major  steps: 
(i)  absorption  of  laser  light,  and 
(ii)  conversion  of  absorbed  optical  energy  to  heat. 

We  consider  strongly  absorbing  semiconductors.  In  strongly  absorbing 
semiconductors,  the  dominant^  absorption  mechanism  at  frequencies  higher 
than  the  bandgap  frequency  is  interband  transitions.  The  conversion  of 
absorbed  energy  to  heat  is  accomplished  through  thermal ization  of 
photoexcited  carriers.  Once  the  conversion  has  been  accomplished,  the 
temoerature  distribution  is  determined  by  solving  the  thermal  diffusion 


equation  in  the  semiconductor.  At  issue  is  the  proper  representation 
of  the  source  term  in  the  heat  equation. 

In  the  model  used  in  this  work,  we  categorize  the  thermal ization 
as  consisting  essentially  of  two  steps,  operative  either  in  parallel  or 
in  sequence:  relaxation  of  the  electron  "excess"  energy  (above  the  bandgap) , 
and  relaxation  of  the  bandgap  portion  of  the  energy.  This  categorization  is  con¬ 
venient  because  the  electron  redistribution  inside  the  conduction  band  proceeds 
generally  at  a  different  rate  than  the  electron-hole  recombination  across  the 
semiconductor  bandgap.  The  bandgap  portion  of  the  electron  energy  is  released 
to  heat  at  the  rate  =  E^  Anr*"^  ,  where  t  is  the  carrier  lifetime,  E^  is 
the  semiconductor  bandgap,  and  An  =  An(x,t)  is  concentration  of  the 
photoexcited  carriers.  An  is  calculated  by  solving  the  electron  diffusion 
equation.  Thus,  the  bandgap  portion  of  the  electron  energy  is  released 
into  a  volume  bounded  by  a  depth  equal  to  max  (L,  a'1),  where  L  *  /Dr  ,  D  is 
carrier  diffusivity  and  a  is  optical  absorption  coefficient. 

The  conversion  of  the  excess  portion  of  the  electron  energy  is 
assumed  in  this  work  to  be  instantaneous.  This  is  because  the  character¬ 
istic  "quasi-thermalization"  frequency  is  on  the  order  of  the  lattice  vibra- 
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tional  frequency,  f  -  10  -  10  sec  ,  which  is  much  faster  than  any 

other  rate  in  the  problem.^ 

The  presence  of  defects  and  interfaces  also  affects  the  heating 
profiles.  Interfaces  affect  heating  mainly  because  of  surface  recom¬ 
bination  of  the  photoexcited  carriers.  The  effect  of  the  surface  will 
be  taken  into  account  in  this  work  through  the  surface  recombination 


Bulk  defects  may  be  categorized  as  localized  or  distributed.  Both 
affect  optical  absorption  as  well  as  carrier  diffusion.  The  influence 
of  the  distributed  defects  will  be  taken  into  account  in  that  the  values 
of  the  optical  absorption  coefficient,  carrier  lifetime,  and  diffusivity 
will  be  different  for  heavily  damaged  than  for  good  quality  materials. 
However,  the  crystal  will  be  considered  homogeneous,  and  its  properties 
independent  of  time.  In  this  sense,  the  values  of  the  material  para¬ 
meters  used  are  some  "average"  values  of  the  parameters  during  the  pro¬ 
cessing.  The  main  reason  for  this  is  mathematical  simplicity.  The 
localized  bulk  defects  will  not  be  treated  here  explicitly.  However, 
their  effect  on  the  heating  should  be  similar  to  that  of  the  surface 
defects.  This  is  because  both  surface  and  localized  defects  are  fast 
recombination  centers,  with  the  subsequent  release  of  the  carrier  energy. 
This  leads  to  the  formation  of  thermal  spikes  at,  and  a  possible  disso¬ 
lution  of,  localized  defects. 

As  has  already  been  mentioned,  the  "bandgap  energy  source  term"  in 
the  heat  equation  requires  knowledge  of  the  electron  concentration,  n(x,t). 
For  constant  T ,  D,  and  v  ,  we  have  solved  the  electron  diffusion  equation 
exactly  for  a  square-shape  laser  pulse  (no  drift  current) .  The  solution 
for  n(x,t)  is  a  closed-form  expression  consisting  of  complementary  error 
functions.  To  the  best  of  our  knowledge,  this  solution  has  never  before 
been  reported  in  the  literature.  The  existence  of  a  closed-form  expression 
for  n(x,t)  is  important  for  our  purposes  since  it  avoids  numerical  inte¬ 
gration  of  the  electron  diffusion  equation.  The  heat  transport  equation 
is  then  solved  numerically,  with  the  temperature  dependence  of  the 
thermal  conductivity  included. 
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In  Section  II  of  this  paper,  we  describe  the  mathematical  basis 
of  the  model.  In  Section  III,  we  discuss  the  results  of  the  numerical 
computations  of  the  temperature.  Section  IV  summarizes  the  results. 


II.  MODEL 

The  photon-to-thermal  energy  conversion  consists  of  photoexcitation 
and  thermal ization  of  the  photoexcited  carriers,  and  heat  transport. 

A.  Concentration  of  Photogenerated  Carriers 

The  instantaneous  concentration  of  the  photogenerated  carriers  is, 
in  the  one-dimensional  approximation,  and  in  the  absence  of  electric 
fields,  obtained  by  solving  the  electron  transport  equation, 


3n 

3t 


M  +  p  i2" 
T  3x^ 


(1) 


where  the  generation  function  g  gives  the  photogeneration  rate  of  the 
carriers , 
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=  aNe 


-ax 


(2) 


where  a  is  the  optical  absorption  coefficient,  and  N  =  N(t)  is  the  photon 
flux  incident  at  the  surface.  N ( t )  =  F(t)/hv,  where  F(t)  is  the  absorbed 
laser  power  and  hv  is  the  photon  energy. 


An  =  n(x,t)  -  nQ  ,  (3) 

where  nQ  is  the  equilibrium  concentration,  x  and  D  are  minority  carrier 
lifetimes  and  diffusion  constant,  respectively. 

In  Equation  (l),a,  r,  and  D  are,  in  general,  functions  of  both  x  and  t. 
The  x-dependence  arises  due  to  the  spatial  inhomogeneities  in  the  semicon¬ 
ductor;  the  time  dependence  arises  from  the  dependence  of  a,  D,  and  x  on 
temperature.  Add i t i  >nal  time  dependence  arises  due  to  changes  in  the 


semiconductor  properties  as  the  damage  disappears  during  the  annealing. 

For  high  carrier  concentrations,  a  may  be  a  function  of  the  photocarrier 
concentration  n(x,t).  It  would  be  quite  complicated  to  fully  account  for 
the  effect  of  the  changing  material  properties.  We  approximate  material 
properties  by  some  "average"  values  and  the  crystal  is  assumed  homogeneous, 
with  the  exception  of  the  front  surface. 

Taking  into  account  surface  recombination  of  the  carriers,  the 
boundary  conditions  are 


and 


j ( x  =  0,t) 


=  -  vs(n(0)  -  n0) 
x=0 


n(x  =  2  )  =  nQ 


(4) 


(5) 


where  v$  is  the  surface  recombination  velocity  and  l  is  the  semiconductor 
thickness.  Equations  (4)  and  (5)  correspond  to  a  finite  flux  of  carriers 
at  the  semiconductor  front  surface,  and  a  constant  carrier  concentration 
at  the  slab  back  boundary.  The  initial  condition  is 


n( x ,  t  =  0)  =  nQ 


(6) 


In  general,  if  N,  a,  t,  and  D  are  functions  of  x  and  t,  Equation  (1) 
has  to  be  solved  numerically.  However,  Equation  (1)  can  be  solved 
exactly  analytically  for  the  special  case  when  a,  t,  and  D  are  constants, 
and  the  generation  function  N(t)  is  a  square-shaped  laser  pulse, 

N(t)  =  NQ  0  <  t  <  t 


N(t)  =  0 


otherwise 


(7) 


In  Eq.  (5),  t  is  the  length  of  the  laser  pulse  and  NQ  is  the  number 
(constant)  of  arriving  photons  at  the  semiconductor  front  surface.  The 
mathematical  methods  are  discussed  in  Appendix  I.  The  result  is 


where  tQ  =  max  (0,  t  -  tp) ,  the  diffusion  length  L  =  /Dt,  and  erf(z)  and 

erfc(z)  are  error  and  complementary  error  functions,  respectively. 

Equation  (8)  is  evaluated  between  the  indicated  upper  and  lower  boun¬ 
daries,  t  and  t  ,  respectively. 

The  importance  of  the  existence  of  solution  (8)  for  our  purposes 
derives  from  eliminating  the  need  for  numerical  integration  of  Eq.  (1). 
n(x,t)  appears  as  a  source  term  in  the  thermal  dif*  'on  equation  dis- 
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cussed  in  the  next  section. 


B.  Thermal  Transport 


The  temperature  distribution  in  the  semiconductor  after  the  applica¬ 
tion  of  the  laser  pulse  is  obtained  by  solving  the  thermal  equations, 

U  _  3_  4  an  =  S(x,t)  # qx 

3t  3X  \C  3 x/  C  ’  '  ' 

\  P  /  P 

where  Cp  and  <  are  the  specific  heat  and  thermal  conductivity,  respectively, 
and  where  the  source  function  S  is  the  heat  flux  deposited  in  the  semicon¬ 
ductor  lattice  per  unit  length.  Equation  (9)  is  solved  subject  to  appro¬ 
priate  boundary  conditions.  In  many  instances,  these  are  (i)  initial  tem¬ 
perature,  (ii)  thermal  flux  at  the  front  surface,  and  (iii)  thermal  flux 
or  temperature  of  the  heat  sink  at  the  back  of  the  sample. 

(!)  Lightly  Damaged  Materials 

Following  our  previous  discussion,  we  express  S  as  a  summation  of 
three  terms,  ,  S2  and  S3.  is  the  "excess"  portion  of  the  electron 
energy,  S2  is  the  bandgap  portion  of  the  electron  energy,  and  S3  is 
energy  released  by  carriers  recombined  at  the  semiconductor  surface. 

The  excess  portion  of  the  electron  energy  is  released  during  quasi- 
thermal ization  of  photoexcited  electrons  inside  the  conduction  band.  Since 
the  characteristic  frequency,  1012  -  1013  sec-1  is  much  faster7  than  any 
other  rate  in  the  problem,  we  approximate 

S-, (x ,t)  =  (hv  -  Eg)a  N( t)e~aX  ,  (10) 

where  v  is  the  laser  frequency  and  Eg  is  the  bandgap  energy. 

The  bandgap  portion  of  the  electron  energy  is  transferred  to  heat  at 
rate  iQ  \  where  tq  is  the  nonradiative  lifetime  of  the  electron. 
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S3(x,t)  =  Egvs(n(x,t)  -  nQ)6(x)  ,  (12) 

where  6(x)  is  the  delta  function.  Mathematically,  the  effect  of  S3(x,t) 

will  be  to  introduce  a  finite  derivative  aT/ax  at  x  =  0;  see  Eq.  (14) 
in  what  follows. 

Similarly  as  for  the  electron  diffusion,  we  consider  an  essen¬ 
tially  uniform  crystal.  It  is  worth  noticing,  however,  that  the 
effect  of  the  presence  of  localized  carrier  "sinks"  in  the  crystal  on 
temperature  could  be  included  through  source  terms  of  the  type  of  Eq.  (12) 
(tfhich  in  turn  introduce  discontinuities  in  3T/ax  ).  Physically,  this 
corresponds  to  a  localized  deposition  of  the  energy,  and  a  subsequent 
formation  of  thermal  spikes  at  defect  sites.  This  makes  dissolution  of 
localized  defects  feasible.  In  this  sense,  location  x  =  0  of  this  paper 
may  be  viewed  as  a  defect  site,  with  the  recombination  rate  for  the 
carriers  equal  to  vsAn(x,t)  5(x). 

Neglecting  loss  of  heat  irradiated  from  the  front  surface  (it  has 
been  estimated  to  be  typically  six  orders  of  magnitude  smaller  than  the 
incoming  laser  flux),  Eq.  (9)  is  then  solved  subject  to  the  following 
initial  and  boundary  conditions: 

T(x,t  *  0)  =  T0  (13) 

?T-(-x--=  0,tl  =  -  la  v  (n(0,t)  -  n  )  (14) 

ax  <  s  o 


T(x  =  l,  t)  =  Tq  (15) 

where  To  is  the  temperature  of  theheat  sink  at  the  back  of  the  sample.  Equation 
(14)  was  obtained  by  substituting  Eq.  (12)  into  Eq.  (9)  and  integrating  the 
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resulting  equation  across  the  front  boundary.  Equation  (9)  has  to  be 
solved  numerically. 


(2)  Heavily  Damaged  Materials  and  K,cp  =  Constants 

For  heavily  damaged  materials  and  K,cp  =  constants,  Eq.  (9)  has  a 
closed-form  solution.  This  is  because  for  extremely  short  carrier  lifetimes, 
rvs  a  <  1  and  TDa  <  1  (i.e.,  bulk  recombination  is  faster  than  either  the 
surface  recombination  or  carrier  diffusion),  the  source  functions  simplify  to 

(n(x  ,t)  -  n  ) 

Vx’t)  5  Eg  -  *  Eg  g  =  EgaF(t)  e"aX  .  (16) 

and 

S3(x,t)  =  0  •  (17) 

This  yields  the  source  term 

S  s  S1  +  S2  =  hv  aN(t)  e_clX  (18) 

and  the  boundary  condition 


3j.(x  =  O.t)  _  n 

3x  u 


(19) 


The  solution  to  Eq.  (9)  is  then  obtained  by  taking  the  limit 
1/t  0  in  Eq.  (A. 10)  of  Appendix  (A. 2).  This  yields 


'J(x’t)  =  7^ 


<  c.  1 


E 

(+.-) 


e±aX  erfc 


N?* 


/(4</cp)"  t' 


(20) 


exp 


(4</cp)f 


+  2r*x  erf( 


TOT 


where  Fq  h\  Nq,  and  the  upper  and  lower  boundaries  are  t' 


=  t 


and 


tQ  =  max(0,t  -  t  ),  respectively. 
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III.  NUMERICAL  RESULTS 


Equation  (9)  was  solved  numerically  for  ruby  and  Nd:YAG  lasers, 

GaAs  and  silicon,  and  various  values  of  laser  and  material  parameters. 

The  values  of  the  material  parameters  ranged  from  those  appropriate 

for  heavily  damaged  materials  to  those  corresponding  to  high-quality 

3 

materials.  The  specific  heat  was  taken  as  constant,  c  =  1.6  J/cm  for  Si 
2 

and  cp  -  1.78  J/cm  for  GaAs.  The  thermal  conductivity  k(T)  was  repre¬ 
sented  by  a  summation  of  exponentials  chosen  to  fit  data  measured  for 

g 

pure  GaAs  and  silicon.  High  impurity  concentrations  tend  to  lower 

these  values;  this  was  not  taken  into  account.  The  carrier  lifetimes 

were  taken  t  =  t  . 

o 

The  laser  pulse  lengths  were  taken  to  be  10  nsec  for  ruby,  and 
100  nsec  for  Nd:YAG  laser  (these  are  the  optimum  lengths  for  heavily 
damaged  materials).  For  GaAs,  the  optical  absorption  coefficient 
(ruby  only)  was  a  =  2  x  104  cm’1.  The  optical  absorption  coeffi¬ 
cient  for  ruby  irradiation  was  ot  =  2  x  10^  cm’1  and  2  x  104  cm’1  for 
crystalline  and  amorphous  silicon,  respectively.  For  Nd:YAG  irradiation, 
n  was  taken  as  a  =  10J  cm  for  amorphous  and  heavily  damaged,  0  ~  500  cm 
for  lightly  damaged,  and  a  =  100  cm’1  for  good  quality  silicon.”  The  values 

of  D  and  T  ranged  from  30  cm^  s’1  to  3  cm^  s’1,  and  10~5  sec  to  10’1®  sec, 
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respectively,  for  silicon,  and  from  200  cm  s  to  15  cm  s’  and  10  sec 
to  10’11  sec,  respectively,  for  GaAs.  The  surface  recombination  velocity 
was  v$  =  2  x  107  cm  s"1 ,  2  x  104  cm  s’1,  and  0. 

Samples  of  the  numerical  results  are  shown  in  Figures  1-12.  The 
profiles  are  shown  as  a  function  of  the  normalized  distance  «x  away  from 
the  semiconductor  front  surface,  with  time  (in  fractions  and  multiples  of 
the  pulse  length  t  )  as  a  parameter.  The  symbols  shown  (in  cgs  units 


unless  specified  otherwise)  have  their  usual  meaning.  The  normalized 


temperature  (T-T0)/(F0/K0a)  is  dimensionless,  is  a  quantity  (equal  to 
one)  which  has  the  dimension  of  thermal  conductivity.  FQ  is  the  actual 

o 

absorbed  power  (in  W/cm  ) .  Most  of  the  profiles  shown  correspond  to 
ATmax  =  ^00  ’  ^ >400°C  (the  exact  number  can  be  calculated  from  para¬ 
meters  specified  in  each  figure).  Before  inspecting  the  profiles,  it  is 
convenient  to  qualitatively  summarize  the  main  results  of  these  calcula¬ 
tions.  This  will  aid  understanding  when  inspecting  the  profiles. 

A.  Heavily  Damaged  Materials 

For  heavily  damaged  materials,  the  effects  of  carrier  diffusion  and 
surface  recombination  on  heat  generation  can  be  neglected.  The  simulated 

profiles  approached  those  calculated  using  Eq.  (20)  (Eq.  (20)  assumes 
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k  =  constant,  rvsD  <  1  and  rDa  <  1).  As  compared  to  Eq.  (20),  the 
effect  of  k  -  <(T)  on  the  profiles  was  to  slightly  increase  the  surface 
temperature,  and  narrow  the  profile  width.  This  is  because  k(T)  decreases 
with  temperature  increasing. 

Surface  Temperature  and  Profile  Width;  For  the  pulse  length  shorter  than 

the  thermal  diffusion  time,  tp  <  tth,  where  tth  =  cp/ica  ,  the  profiles  are 

optical  absorption  depth-limited.  The  profile  width  «  cT^ ,  and  aT  grows 

max 

linearly  with  both  the  pulse  length  and  optical  absorption  coefficient, 

Almax  '  Foy/y  For  the  Pu1se  length  longer  than  tth,  tp  >  tth,  the 
profile  is  thermal  diffusion-limited,  its  width  is  =  («<:tp/cp) 1  ,  and 

ATmax  is  independent  of  a  and  proportional  to  the  square  root  of  the  pulse 
length,  i’nlax=  Fq  /tp/<cp  (the  surface  temperature  approaches  saturation 
for  extremely  long  pulses).  In  either  case,  ATmgx  is  approximately  linear 
with  pulse  peak  power  (.AT^  is  exactly  linear  with  Fq  for  k  =  constant).  The 
heating  is  most  efficient  (least  laser  energy  required)  when  t  =  t^. 


The  optical  depth-limited  profiles  are  generally  narrower  than  the 
thermal  diffusion-limited  profiles. 

Temperature  Decay:  After  the  end  of  the  laser  pulse,  the  temperature 

drops  and  the  profile  broadens  as  the  heat  propagates  into  the  bulk.  The 

decay  is  much  faster  for  strong  absorption  than  for  weak  absorption.  This 

2 

is  because  the  thermal  diffusion  time  cp/ica  is  much  longer  for  small  a  than 
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for  large  a.  For  example,  for  a  =  2  x  10  cm  ,  the  surface  temperature 
typically  dropped  to  one-half  (1/2)  of  its  maximum  value  at  t  =  100  t  . 

For  a  =  2  x  10^  cm  \  the  halftime  was  only  t  =  5  t  . 

For  comparison  purposes,  Eq.  (9)  was,  for  heavily  damaged  materials, 
also  solved  for  the  "Gaussian"  pulse  shapes.  We  have  observed  really  no 
significant  difference  between  the  shapes  of  the  thermal  profiles  created 
by  square  and  Gaussian  pulse  shapes.  The  resulting  temperature  distribu¬ 
tions  were  affected  much  more  by  the  condition  a”1  vs.  4:tp/cp,  than  by 
the  pulse  shape.  For  the  Gaussian  pulse,  slightly  higher  laser  powers 
were  required  to  achieve  the  same  surface  temperature  as  for  the  square- 
shape  pulse.  The  Gaussian  profiles  tend  to  decay  more  slowly  than  the 
square-shaped  pulse  profiles. 

B.  Lightly  Damaged  Materials 

The  main  effect  of  carrier  diffusion  on  temperature  profiles  consists 
of: 

(i)  Higher  laser  powers  are  required  to  raise  the  surface  tempera¬ 
ture  of  a  good-quality  material  than  of  a  damaged  material; 

(ii)  Temperature  profiles  in  good  quality  materials  are  generally 
wider  than  in  damaged  materials; 

(iii)  The  time  (t  )  at  which  the  surface  temperature  reaches 
maximum  is  longer  in  good  quality  materials  than  in  damaged  materials; 
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(v)  The  effect  of  carrier  diffusion  on  temperature  profiles  is 
much  stronger  for  Nd:YAG  laser  than  for  ruby  laser  irradiation. 

( vi )  The  effect  of  the  linking"  of  the  carriers  is  to  produce 

thermal  spikes  at  the  surface  and  lower  temperature  in  the  material  interior. 
For  a  sufficiently  fast  surface  recombination,  the  profiles  may  develop  a 

secondary  local  maximum  in  the  material  interior. 

The  conclusions  (i)  -  (v)  are  the  direct  consequence  of  the  model 

used  for  the  transfer  of  the  photon  energy  to  heat.  The  excess  portion  of 

the  electron  energy  is  transferred  during  the  "quasi-thermal i- 

zation"  of  the  photoexcited  carriers,  into  the  volume  bounded  by  a"1.  The 

corresponding  heated  volume  is  o  ,, 

eff 

°eff  =  max  («_1 , (< tp/cp)1/2)  ,  (21) 

and  the  associated  rise  in  the  surface  temperature  is  on  the  order  of 


AT  (x=0, 
max'  ’ 


t=V 


.  Vp  . 

Cp°eff 


(22) 


» 


The  bandgap  portion  of  the  energy  heats  a  volume  bounded  by  the  depth 


i 

eff  ,g 


=  max(oeff,L) 


(23) 


where  L,  L  =  /~Dt ,  is  the  carrier  diffusion  length, 
ture  rise  is  roughly 


The  associated  tempera- 


ATmax  (x=0’  t=V  "  c_£ 


F°y 


P  eff,g 


ia 

hv 


I 
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Thus,  in  good  quality  materials  where  L  >  °eff>  the  bandgap 
energy  heats  a  volume  which  is  larger  than  the  one  heated  by  the  excess 
energy.  This  results  in  temperatures  which  are  generally  lower  (for 
the  same  laser  power)  than  those  in  damaged  materials.  Furthermore,  since 
both  the  electron  lifetime  (t)  and  the  thermal  diffusion  time  associated 
with  the  bandgap  portion  of  the  energy 


tth,g 


cp  *eff,g/4* 


(25) 


may  be  quite  long;  the  heating  times  (tffl  and  t-j  ^  may  be  substantially 
longer  than  tp.  Clearly,  the  effect  of  carrier  diffusion  on  the  tempera¬ 
ture  will  be  greater  for  NdrYAG  than  for  ruby  laser  irradiation.  This  is 
because  Nd:YAG  frequency  is  very  close  to  the  silicon  bandgap  frequency;  the 
excess  portion  of  the  electron  energy  is  =  1%.  For  ruby  irradiation  of 
silicon,  the  excess  portion  of  the  electron  energy  is  large,  =  35%,  and 
is  not  affected  by  electron  diffusion.10 


The  conclusion  (vi)  is  the  unique  consequence  of  the  surface  recom¬ 
bination  of  the  carriers.  For  v$  =  0,  the  temperature  is  smooth  and 
monotonically  decreasing  with  the  distance  away  from  the  surface.  This 
is  because  the  thermal  sources  (optical  absorption)  and  S£  (carrier  dis¬ 
tribution)  are  also  monotonically  decreasing,and  the  surface  recombination  source 
=  0.  For  vs  f  0,  however,  the  carriers  preferentially  recombine,  and  release 
their  energy,  at  the  surface.  This  leads  to  the  formation  of  thermal  spikes. 
Furthermore,  since  the  carrier  distribution  develops  a  "sink"  at  the  surface, 
the  source  S2(x)  will  be  nonmonotonic.  If  the  maximum  is  sufficiently 
strong  so  that  also  the  total  S  =  S-j  +  is  nonmonotonic,  then  the 
resulting  thermal  distribution  may  develop  a  local  interior  maximum. 
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(We  note  that  all  peaks  and  maxima  will  eventually  disappear  with  time 
because  of  carrier  and/or  thermal  diffusion.) 


Figures  1  and  2  show  carrier  distributions  for  silicon,  with  time 
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(t/tp)  and  v$  as  parameters.  In  Figure  1,  the  electron  lifetime  is  10" 

sec,  in  Figure  2,  t  =  10"7  sec.  The  effect  of  v$  is  clearly  seen.  The 

carrier  distribution  (per  unit  photon  flux)  is  high  and  monotonic  for 

vs  =  0,  lower  and  nonmonotonic  for  v$  f  0.  For  short  t  (Figure  1),  the 

maximum  is  sharp  and  close  to  the  surface.  The  maximum  moves  into  the 

material  interior  with  increasing  v$  and  higher  t/tp.  Because  of  the 

short  t  (t  <  tp) ,  the  distribution  decays  extremely  rapidly  after  the 

end  of  the  pulse  (t  >  tp) .  For  longer  T (t  l  tp,  see  Figure  2),  the 

maximum  is  much  more  pronounced.  The  concentration  remains  high  even 

when  t  >  t  . 

P 

Inspection  of  Figures  3-12  reveals  the  qualitative  features  discussed 
above. 

Ruby  Irradiation  of  GaAs 

Figures  3,  4,  and  5  show  temperature  distributions  for  ruby  laser 
irradiation  of  GaAs.  Figures  3  and  4  show  the  effect  of  v$.  The  principal 
effect  of  large  v$  in  GaAs  is  to  lower  the  surface  temperature.  For  the 
same  laser  input,  the  surface  temperature  obtained  for  vg  =  2  x  104  cm  s'1 

is  higher  (by  =  50%)  than  the  temperature  obtained  for  v$  =  2  x  107  cm  s“^ . 

Figures  4  and  5  compare  temperature  distributions  for  heavily  damaged  and 
good  quality  GaAs.  The  principal  effect  of  the  electron  diffusion  is  to 
increase  the  laser  power  required  for  heating  good  quality  material  in 
Figure  5,  approximately  by  (1  -  Eg/hv)’^  -  5  (for  the  same  a).  Its  effect 


on  the  width  of  the  profiles  is  negligible;  the  widths  of  profiles  shown 
in  Figures  4  and  5  are  practically  identical.  This  is  because  in  good 
quality  material,  the  bandgap  energy  (80%  of  the  photon  energy)  is 

O 

distributed  into  the  depth  L  -  3  x  10  cm,  which  is  much  larger  (by 
nearly  two  orders  of  magnitude)  than  the  depth  for  the  excess  energy 
(20%),  a  -  5  x  10  cm.  This  makes  the  bandgap  portion  of  the  energy 
ineffective  for  the  heating.  We  should  note  that  for  intermediately 
damaged  GaAs ,  more  of  the  available  power  will  be  used  for  the  heating, 
and  for  -  L,  the  profiles  will  be  somewhat  wider  than  in  either  the 
heavily  damaged  or  high  quality  materials. 

Ruby  Irradiation  of  Silicon 

Figures  6-7  show,  respectively,  temperature  distributions  in  lightly 
damaged  and  good  quality  silicon  irradiated  by  ruby  laser  radiation.  We 
have  also  simulated  profiles  in  amorphous  silicon.  The  power  requirements 

and  shapes  of  these  profiles  were  very  much  the  same  as  those  shown  for 
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GaAs.  This  was  because  of  the  similar  input  parameters  (ct  =  2  x  10  cm  , 

_g 

t  5  10  sec)  used  in  these  calculations. 

Figure  6  corresponds  to  a  lightly  damaged  crystalline  silicon 
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(a  =  2  x  10  cm"  ).  The  main  di fferences between  amorphous  and  crystalline 
silicon  are  higher  power  requirements  and  longer  heating  times  for 

crystalline  materials.  The  higher  powers  are  required  because  of 

(i)  weaker  optical  absorption,  and  (ii)  inefficient  utilization  of  the 

bandgap  portion  of  the  energy  for  the  heating  in  good  quality  materials. 

The  heating  times  (t^2  =  100  t  in  Figure  6)  are  long  because  of  the 
influence  of  the  bandgap  energy  on  temperature  (for  T  >  tp)  at  long 


times  (in  Figure  6,  t  =  10  t  ) .  At  short  times,  the  profile  width 
is  -  a"\  as  only  the  excess  energy  heats  the  material  at  t  <  t. 

Figure  7  corresponds  to  a  good  quality  material  (t  =  10'^  sec, 
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D  =  30  cm  s"  ).  The  interesting  feature  to  notice  in  Figure  7  is 
that  at  times  t  =  100  tp,  the  profile  width  is  narrower  than  the  one 
shown  in  Figure  6.  This  is  because  at  long  T  and  D  the  bandgap  energy 
has  practically  no  effect  on  the  heating  (at  sufficiently  short  times). 

For  L  >  cf \  both  ATmax  and  the  halfwidth  are  governed  by  the  excess 
portion  of  the  electron  energy,  the  bandgap  portion  affects  the 
profile  only  at  very  long  times,  t  *  t. 

Nd:YAG  Irradiation  of  Silicon 

Several  representative  profiles  are  shown  in  Figures  8-12.  For 
Nd:YAG  in  silicon,  Eg  -  .99  hv  and  the  profiles  are  therefore  bandgap 
recombination  dominated.  Figures  8  and  9  compare  profiles  for  two 
values  of  v$  for  heavily  damaged  material  (a  =  103  cnf\T  =  10~9  sec); 
vs  =  2  x  107  cm  s"1  and  2  x  104  cm  s"1  in  Figures  8  and  9,  respectively. 

The  halfwidths  are  =  >  L  in  both  cases.  The  primary  effect  of  slow 

v  (v  t  < L)  in  Figure  9  is  to  slightly  increase  AT  and  suppress  the 
interior  maximum  appearing  in  Figure  8  for  large  v$.  There  should  be 
a  sharp  peak  at  the  surface  in  Figure  8.  The  peak  is  unresolved  because 
of  the  numerical  limitations.  The  peak  corresponds  to  thermal  spikes 
(discussed  earlier),  created  at  x  =  0  as  the  result  of  the  preferential 
deposition  of  the  carrier  energy  at  the  surface.  Several  remarks  about  the 
effect  of  vs  f  0  on  the  profiles  are  in  order  here. 

As  has  already  been  discussed,  the  occurrence  of  thermal  spikes 
and  interior  maxima  (the  local  maximum  is  always  lower  than  the  spike) 


is  the  consequence  of  fast  surface  recombination.  By  inspecting 
the  simulated  profiles,  a  nonmonotonic  T(x)  occurs  for:  v$  suffi¬ 
ciently  large,  vst  >  L  (so  that  maximum  in  n(x,t)  is  large  and  separated 
from  the  surface),  t  not  too  short  (so  that  the  energy  does  not  relax 
instantaneously),  a_1  not  too  small,  a_1  £  L  (so  that  the  bandgap  energy 
is  not  dissipated  into  too  large  a  volume),  and  the  excess  portion  of  the 
energy  not  too  large  compared  to  the  bandgap  energy  (so  that  the  source 
term  S^(x)  does  not  suppress  the  peak  in  S£(x)).  If  these  conditions  are 
satisfied,  T(x)  becomes  nonmonotonic,  with  a  sharp  spike  at  the  surface 

surface  and  a  local  maximum  at  the  distance  0  <  x  <  L  away  from 
the  surface.  We  note  that  as  the  carriers  diffuse  deeper 
into  the  material,  the  local  maximum  may  also  move  deeper  into  the 
interior;  it  eventually  disappears  (  because  of  both  decreased 

carrier  concentration  and  thermal  diffusion),  with  the  characteristic 

2 

time  t  =  max  (T,  t_,  (c  /4<)L  .  The  magnitude  of  the  spike  at  x  =  0 
P  P 

can  be  estimated  by  using  the  energy  conservation  arguments.  For  example, 
the  energy  released  by  the  carriers  recombined  at  the  surface  must  be 
equal  to  the  heat  contained  in  the  spike.  We  then  obtain  that  in  the 

limit  of  rapid  vs ,  short  t,  and  large  a,  the  height  of  the  spike  is  on 

the  order  of  AT  .  =  aiF  ,  where  I  is  the  characteristic 

peak  o  s  s  p  s 

1  ength associated  with  the  position  of  the  maximum  in  An(x,t).  For 

O 

typical  parameters  of  interest,  this  yields  ATpgak  =  several  10  °K  and 

spike  characteristic  width  z  =  /?.  Kfc  v  =  10'6  cm. 

p  s  p  s 

Figure  10  corresponds  to  light-to-intermediately  damaged 
material  (a  *  103  cm"1,  t  -  10'7  sec).  The  development  of 
nonmonotonic  AT(x,t)  is  clearly  seen  (again,  there  should  be  a  spike 

at  x  =  0).  The  temperature  becomes  monotonic  at  times  t  X 
7  6 

ttu  ~  (c  /4*)L  -  10  sec.  Because  of  the  long  t,  the  required 

th,e  p 


power  is  about  twice  higher  than  the  one  for  the  heavily  damaged  material 

in  Figure  9  (for  otherwise  the  same  parameters).  Figures  11  and  12  show 

the  effect  of  the  electron  diffusivity  D  on  the  profiles.  The  main 

effect  of  increasing  D  (at  constant  r,  t  ~  10"^  sec)  is  to  suppress  the 

interior  maximum.  Because  of  the  long  t,  the  heating  times  are  very  long, 

the  temperature  is  high  even  at  time  t  £  100  tp. 

We  have  also  simulated  profiles  for  a  =  100  cm"^  ,  and  500  cm”^ .  On 

the  scaled  distances  ax,  the  shapes  of  the  profiles  were  qualitatively  similar 
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to  those  calculated  for  a  =  10  cm  .  The  required  powers  were  also  much 

higher  (approximately  by  the  ratio  of  a's),  and  the  heating  times  accord¬ 
ingly  longer  (as  determined  by  the  thermal  time  (cp/4*)j^ff  see  Eq.  25). 

IV.  DISCUSSION 

In  order  to  facilitate  a  more  comprehensive  understanding  of  our 
results,  Figures  13  and  14  show  the  derived  dependence  of  (i)  pulse 
energy  (Q)  required  to  raise  the  surface  temperature  by  aT  =  1  ,400°C, 

(if)  time  at  which  the  surface  temperature  reaches  its  maximum  (tm),  and 
(iii)  time  at  which  the  temperature  drops  to  one-half  of  its  maximum 
value  (t]/2),  on  the  quality  of  the  processed  material  (silicon  only). 

The  quality  is  characterized  in  terms  of  the  electron  diffusion  length  L 
(the  corresponding  values  of  x  and  D  and  a  are  also  indicated).  We 
emphasize  once  again  that  the  values  of  the  material  parameters  used  are 
some  "average"  representative  values  (this  is  true  particularly  for 
the  values  of  a  for  Nd :YAG  irradiation  of  silicon)  and  thus  the 
values  of  the  derived  parameters  are  very  approximate.  The  derived 


QlQ  =  F  tp)  Is  the  actually  absorbed  energy:  the  Incident  energy  is 
nearly  twice  as  high  as  the  absorbed  energy,  because  of  material  reflection. 
The  laser  parameters  have  the  same  values  as  those  used  in  Section  III; 
the  results  could  be  generalized  to  different  values. 

A.  Ruby  in  Silicon 

The  required  energy  increases  with  L  increasing,  from  about 
Q  -  .25J  for  amorphous,  to  Q  =  4  J  for  good  quality  material.  Only  the 
excess  portion  of  the  energy  is  operative  for  heating  of  the  quality 
material  (Eg/hv  -  65%),  aind,  therefore,  Q  saturates  at  long  L;  this  is 
seen  in  Figure  13.  With  the  a's  and  tp  used,  the  heated  volume  is 
thermal  diffusion-dominated  for  amorphous  material  and  optical  absorption 
depth-limited  for  good  quality  material.  The  surface  temperature  reaches 
maximum  at  tm  =  tp.  The  effect  of  the  bandgap  portion  of  the  energy  is 
to  slow  down  the  temperature  decay  in  lightly  damaged  materials.  At  Inter¬ 
mediate  values  of  L,  L  *  of1,  the  temperature  remains  high  up  to  the  time 

(cp/4<)L  (i.e.,  t1/2  „  t).  At  very  long  L  (L  »  a-1),  however,  the 

bandgap  energy  is  dissipated  into  an  extremely  large  volume  (depth  = 
several  L),  so  that  its  effect  on  temperature  becomes  insignificant. 

The  temperature  starts  to  decay  more  rapidly.  The  decrease  in  t^2  is 
seen  in  Figure  14  for  L  K  5  x  10"3  cm. 

B.  Nd;YAG  in  Silicon 

The  heating  is  nearly  entirely  bandgap  energy-dominated.  For  amorphous 
and  damaged  materials,  the  optical  absorption  depth  a-1  >  (ath,l),and  the 
profile  width  is  therefore  «  a'1,  unaffected  by  the  diffusion.  The  surface 
temperature  is  ATmax  -  Qa/cp.  In  good  quality  material,  L  becomes  K  10"2  cm 


This  is  larger  than  the  probable  value  of  a  in  the  material, 
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a-1  =  2  x  10-3  cm  for  lightly  damaged,  and  a  -  10" 

cm  for  good  quality  material.  Therefore,  the  depth  of  the  heated  volume 
is  -  L,  and  the  power  required  to  heat  the  material  increases  accordingly. 
In  Figure  13,  Q  increases  at  long  L  with  L  approximately  linearly. 

Because  the  heating  is  bandgap  energy-dominated,  the  maximum  value 
of  the  temperature  in  good  quality  materials  is  obtained  at  time  tm  , 
generally  longer  than  tp  (for  the  amorphous  and  heavily  damaged  material, 
tm  =  tp).  For  good  quality  material,  tm  is  on  the  order  of  several  t. 

The  decay  of  the  surface  temperature  is  determined  by  the  thermal  diffu¬ 
sion  time  tth>g  (t]/2  tth>g  =  max[(cp/4*c)L2,  (cp/4*)  t  ].  Since 
the  thickness  of  the  heated  region  increases  with  L  increasing,  there 
is  no  saturation  of  t^2  at  long  L  which  has  been  observed  for  ruby 
irradiation  of  silicon. 

In  summary,  we  have  investigated  in  this  work  the  influence  of  the 
diffusion  of  the  photoexcited  carriers  on  laser-heating  of  strongly 
absorbing  semiconductors.  The  diffusion  was  found  to  increase  the 
power  required  to  heat  good  qual i ty  materials ,  and  to  prolong  the 
heating  times.  Surface  recombination  was  found  to  create  thermal  spikes 
at  the  surface.  The  temperature  profiles  calculated  in  this  work  indi¬ 
cate  that,  under  certain  conditions,  the  temperature  distribution  may 
become  nonmonotonic.  The  nonmonotonic  profile  is  the  consequence  of 
the  interplay  of  surface  recombination  and  diffusion  of  the  photoexcited 
carriers.  The  position  of  the  interior  maximum  is  dependent  upon  both 
the  bulk  quality  and  the  properties  of  the  surface.  The  strength  of  the 
thermal  spike  can  be  estimated  from  energy  arguments. 

In  practice,  the  existence  of  the  thermal  spikes  may  be  of  importance 
for  annealing  of  localized  defects.  For  a  sufficiently  strong  spike,  it 
is  qualitatively  feasible  for  the  heat  to  dissolve  the  defect.  This 


results  in  localized  defect  removal.  For  practical  utilization, 
further  analyses  are  desirable.  This  is  because  of  the  effect  of  both 
the  finite  size  (e.g.,  finite  thickness  of  the  surface  "dead  layer"),  and 
spatial  distribution,  of  the  defects  on  the  temperature  profiles. 

Further  effects,  not  considered  here,  include  other  energy  "diffusion" 
mechanisms  (such  as  radiative  recombination  of  the  carriers),  non¬ 
equilibrium  effects,  and  distributed  damage  such  as  that  created  during 
ion  implantation.  The  damaged  regions  have  different  optical,  thermal, 
and  electronic  properties  than  the  rest  of  the  material  (they  often 
aabsorb  the  radiation  more  strongly).  Furthermore,  as  the  damage  dis¬ 
appears,  material  properties  change  during  the  processing.  This  was  not 
taken  into  account  here;  in  the  model  used,  with  the  exception  of  the 
front  surface,  the  material  was  assumed  homogeneous,  and  its  properties 
independent  of  time.  Further  analyses  are  being  planned  to  account 
for  these  effects. 
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APPENDIX  A:  ELECTRON  TRANSPORT  EQUATION:  SOLUTION  BY  LAPLACE 
TRANSFORM  TECHNIQUE 


A.  1.  Problem 


The  diffusion  equation  we  want  to  solve  is 

(n  -  O 


3n  c/+\  -cix 
at  aF(t)e 


+  D 


32n 

3x2 


(A.l ) 


We  consider  constant  a,  t  and  D,  and  a  square-shape  generation  function 


F(t): 


0  <  t  <  tp 
otherwise. 


(A. 2) 


The  boundary  conditions  are: 

Front  surface: 

0|V  -  vs(n(0,t)  -  n0)  (A. 3) 

/x=0 

Back  surface: 

n(x  =  l)  =  nQ  (A. 4) 

The  initial  condition  is 


n(x,t=0)  =  nQ  (A. 5) 

Equation  (A.l)  is  a  linear  partial  differential  equation  for  n(x,t). 

For  vs  =  0,  Eq.  (A.l)  can  be  solved  by  using  the  Green's  function  technique. 
For  v$  f  0,  it  is  solved  by  using  the  Laplace  transform  technique.  In  using 
the  Laplace  transforms,  advantage  is  taken  of  the  knowledge  of  the  steady 
state  solution  to  (A.l). 
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2.  Case  v  =0:  Green's  Function  Technique 
s  _ 


For  v$  =  0,  Green's  function  .is 


G(x,x'  ,t-t' )  = 


+oo 

~—r  5>Un 


✓4irD(t-f) 

l ( x-x '  -  2ns.)2\  .  ..._/ (x  +  x'  -  2ns02\\ 

exP  ^  4D(t  I  V)  )  +  6XP  \  4D(t  -  H  /  ]  ’ 


and  An  =  n(x,t)  -  nQ  is 


(A. 6) 


oo  p 

An(x  ,t)  =  J  dx'  J  G(x,x',t-t')  aF0  exp[-ax] 


(A. 7) 


Assuming  the  slab  thickness  l  »  a~\  the  semiconductor  can  be 

replaced  by  a  semi-infinite  medium,  with  the  result  that  the  summation 

l  in  Eq.  (A. 6)  reduces  to  only  one  term,  namely  the  term  n=0.  Sub- 
n 

stituting  the  resulting  Equation  into  Eq.  (A. 7)  and  integrating  it  over 
coordinate  x  yields: 


An(x,t) 


(A. 8) 


exp 


[- 


+  2aDt' 

2/Dtr 


where  L  =  /ih,  erf  c(z)  is  the  complementary  error  function,  and  the 


integral  lower  boundary  tQ  is 


t>  tp 


(A. 9) 


Integral  (A. 8)  is  evaluated  by  using  the  partial  fraction.  This 


yields 


,t)  =  - 2"1 - X]  exp{;ax)  expt^  (l-a2l2)l  erfcfr—  2aClt  -1 

2 ( ot2L2- 1 )  LT  J  L  2/Dt7-  Jy 


'i  /  C-aV)] 


exp(±ax) 


4/D 


( 


The  second  term  on  the  right  hand  side  of  Eq.  (A. 10)  is  evaluated  by 


using  the  following  change  of  variables. 


-?  2 
t  =  z  ,  and  t  =  z 


(A. 11) 


for  the  first  and  second  term  under  the  integral  sign,  respectively.  Per¬ 
forming  the  change  of  the  variables,  each  of  the  two  terms  transforms  into 


an  expression  of  the  type  (a  f  0) 


/ «p(-a2z2-^j  dz  =  §  ^ 


exp±2ab)  erf(az  ±  b)+  constant  .  (A. 12) 


Collecting  the  resulting  terms  completes  the  calculation.  The  final  result 


i(x,t)  =  - -  exp|-{  1  -L2a2)—  j  X!  e±aX  erfc(a/Dt'  ± 

2(a2L2-l)  L  TJ("> 


+  aL  Xi  e  ^  erf  f  - 

(  +  ,-) 


»w 


mr 


(A. 13) 
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A. 3.  Case  v$  f  0:  Laplace  Transform  Technique 

Since  the  steady-state  solution  to  Eq.  (A.l)  is  known,  Eq.  (A.l)  is 
Laplace- transformed  with  respect  to  time.  One  obtains  (for  the  square¬ 
shaped  F( t) ) : 


0  e-*x  (1  -  e'S  p)  -  (s«l)  An(x,s  )  ♦  D  LsLlfL  -  0  , 


(A. 14) 


where  n  =  n(s,t)  is  the  Laplace  transform  of  n(x,t).  The  solution  of 
(A. 13)  is 


C*F  T  n  -5Lp. 

An( x ,  s )  =  - - 0 e  P)  o 

(1  +  St  -  oA2) 


V$  +  aD 


(A-15) 


The  reverse  Laplace  transform  Is  easiest  to  find  by  using  the 
well-known  results  for  the  partial  expressions  appearing  in  Eq.  (A. 15), 

and  making  time  convolutions  of  the  partial  Laplace  transforms.  That  is. 


F(t)  =  f Fj (t  -  t')  F2(t')dt'  , 


(A. 16) 


where 


fi  (s) 


F,(t) 


f2(s) 


^(s)  fz(s) 
f(s  +  c) 


f  (as) 


F2(t) 


F(t)  e' 


The  partial  transforms  used  In  evaluating  the  Inverse  Laplace  transform 
of  Eq.  (A. 15)  are  of  the  type 


1(1  .  e-b(t-a))  H(t_a)  a  >  o 


e~as 

s(s+b) 


(A. 17) 


^0  -  e"bt) 


a  -  0 


and 

-  b  expCab  +  b2t)  erfc^^—  +  b/t^  .  (A. 18) 

In  Eq.  (A. 17),  H ( t)  is  the  step  function.  H(t)  =  1  for  t  >  0  and  zero 
otherwise. 


s±E 

b  +/s 


/irt 


exp 


(■ 


Using  (A. 16)  -  (A. 18),  the  inverse  transform  of  Eq.  (A. 15)  is 


Integral  tA.19)  is  evaluated  by  using  the  same  changes  of  variables  as 
those  used  for  evaluation  of  integrals  in  the  Green's  function  expression 
in  Appendix  A.2(viz.,  Eq.  (A. 8)  and  the  following  discussion).  It  is 
convenient  to  rewrite  Eq.  (A. 19)  in  the  form: 


The  solution  of  Eq.  (A. 20)  is 

An(x 


where 


cxtF„  ( 

;,t) 


exp 


-aX-^-(l-aZL2) 


vs+aD 


x  (F(t*}-(G(t' ) -G(0))exp 


-f-0-<A2) 


F(t)  =  Jit'  f(t')  +  constant 

G(t)  =f dt'  f(t')  g(t')  +  constant 


(A. 2 


(A.  2i 
(A.2‘ 


The  first  term  in  the  integral  on  the  right-hand  side  of  Eq.  (A. 24)  is 
evaluated  by  using  the  change  of  the  variables  t'  =  z  .  The  resulting 
integral  is  of  the  type  of  Eq.  (A. 12).  The  second  term  under  the  integra¬ 
tion  sign  is  evaluated  by  first  using  the  partial  fraction,  and  then  making 

_2 

the  change  of  the  variables,  t'  =  z  in  the  two  resulting  integral 
expressions.  The  integral  is  of  the  type 


/«  =  "b<H>  erfH 1  f)  ■  ( 

Integral  (A. 25)  is  evaluated  by  using  the  same  procedure.  In  fact,  the 
only  difference  between  the  integrals,  Eq.  (A. 24)  and  Eq.  (A. 25),  is  the 
replacement  of  t"  in  Eq.  (A. 24)  by  a  D  in  Eq.  (A. 25). 


(A. 26) 


f(t,t)  g(t,x)  — *■  f(t, — j) 

Da 


(A. 27) 


The  mathematics  is  then  straightforward,  though  rather  tedious.  After 
having  collected  and  combined  all  terms,  the  final  result  is 


«(«.») . 


(1  +  ccD/v  \  ±  f  , 


{±  1  +  D«/\) 

+  /  ^  ~  - -  exp 

,  :  0  -  Da/V  ) 

(+,-)  S 


[±  ax  -  ^l-A2)]  erfc/±  — * —  +  a/Dt\ 

V  mt>  ) 


(A. 28) 


+  _  20  -  L2a2) 


- — 1  —  L"a^ - ~ - exp|^’+  ^  I  erfc(vsJI'  +  ~~ — 

UD/v$  -  1)(1  -  vs2t/D)  [t  V  /*Dt« 


fmax(0,t-tp) 


I 


*1 


It  is  easy  to  verify  that  Eq.  (A. 28)  reduces  to  Eq.  (A. 13)  in  the 
limit  of  v$  =  0.  In  the  steady  state,  Eq.  (A. 28)  reduces  to  the  well-known 
expression , 


Limit  An(x,«) 

tp,  t— 


ctrf  _ 


2.2 
a  L 


vs  +  aD 


-x/L 

e 


(A. 29) 
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9.  Temperature  for  Nd:YAG  irradiation  of  damaged  silicon:  t  =  10'  sec, 
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^  r 
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10. 


Temperature  for  Nd:YAG  Irradiation  of  silicon,  t  =  10"®  sec, 

D  =  3  cm2  s"1 ,  v$  =  2  x  104  era  s"1 ,  a  =  103  cm"1 ,  t  =  10~7  sec, 

FQ  =  2  x  108  W/cm2. 

Temperature  for  Nd:YAG  irradiation  of  silicon,  t  =  10"®  sec, 

D  =  30  cm2  s'1 ,  v  =  2  x  104  cm  s_1 ,  a  =  103  cm"1 ,  t  =  10“7  sec, 

s  p 

Fq  =  2.5  x  108  W/cm2. 

Laser  energy  (Q)  required  to  heat  the  silicon  surface  just  below 
the  melting  temperature.  Nd:YAG:  t  =  10"7  sec;  ruby:  t  =  10"8 

p  J  p 

The  dependence  of  (1)  time  at  which  the  surface  temperature  reaches 
the  maximum  (tm);  time  when  the  surface  temperature  decays  to  1/2 
of  its  maximum  value  (t1/2),  on  material  quality  (silicon). 


01 


the  melting  temperature. 
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Ruby  laser  pulses  of  1.5  to  2.5  J  cm  were  shown  to  be 

effective  in  dissolving  defect  nuclei  in  {100}  and  {111} 

14  15  -2  14 

silicon  surfaces  implanted  with  10  to  10  B  cm  and  10  to 
14  -2 

5  x  10  P  cm  .  No  melting  and  only  partial  electrical  activation 

occurred.  A  sequential  moderate  thermal  anneal,  10  min  at  1000°C, 

results  in  full  electrical  activation  without  significant  impurity 

redistribution.  A  {111}  silicon  surface,  made  amorphous  by  an 

14  -2 

ion  implantation  of  5  x  10  P  cm  was  electrically  activated  and 
completely  cleared  of  defect  nuclei  without  melting  by  preceding 
the  laser  anneal  by  a  two  step  thermal  anneal  (30  min  at  550°C, 

10  min  at  1000°C) .  The  mechanism  of  pulse  laser  annealing  of 
non-amorphous  ion  implanted  surfaces  is  discussed. 


B  2 


LASER  ANNEALING  OF  LOW  FLUENCE 


ION  IMPLANTED  SILICON 


1.  INTRODUCTION 

Ion  implantation  has  been  found  to  result  in  the  intro¬ 
duction  of  electrically  active  impurities  in  an  extremely 
uniform  and  reproducible  manner.  Concomitant  with  this  de¬ 
sirable  characteristic,  ion  implantation  leaves  a  residue 
of  lattice  defects.  These  primary  defects  consist  of  impurity 
atoms  on  non-lattice  positions,  and  the  presence  of  many 
silicon  vacancies  and  interstitials.  A  relatively  moderate 
thermal  anneal,  1000°C  for  10  minutes,  results  in  relocation 
of  the  impurity  atoms  to  lattice  positions  as  shown  by  their 
electrical  activation.  There  is  however  a  considerable 
residuum  of  lattice  damage  remaining.  These  secondary  defects  ^ 
have  been  studied  by  TEM  for  both  light  and  heavy  implantations 

and  have  been  shown  to  consist  primarily  of  dislocation  loops, 

(2 ) 

with  and  without  enclosed  stacking  faults.  For  heavy 

implantations  the  Rutherford  ion-backscattering  and  ion- 
channel  technique  has  been  used  to  show  the  presence  of  residual 
defects.  ^  A  method  developed  by  the  author  ^  utilizes 
a  two  step  technique.  By  growing  a  thermal  oxide,  extrinsic 
dislocation  loops  are  found  to  expand  to  a  size  of  several 
microns.  Using  chemical  etching  techniques  the  inter¬ 
section  of  these  ternary  defects  with  the  surface  can  be 
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observed  by  optical  or  scanning  electron  microscopy.  This 
method  was  applied  to  boron  &  phosphorus  over  the  range  of 
lO"*-2  to  10'1'~’  cm2.  A  thermal  anneal  of  1  hr  at  1100°C  in  dry 
N2  was  effective  for  1014  cm2  but  not  for  1015  cm2  fluences 
in  dissolving  the  defect  nuclei  introduced  by  the  ion  im¬ 
plantation  so  that  a  subsequent  oxidation  failed  to  result 

( 4 ) 

in  the  development  of  expanded  dislocation  loops.  This 

however  causes  a  significant  redistribution  of  the  implanted 
impurities . 

In  this  study  a  Q-switched  ruby  laser  was  used.  A 
review  of  the  literature,  as  well  as  a  recent  review  article 
indicates  that  this  laser  has  only  been  effectively  used  when 
the  implant  fluences  were  high  enough  to  convert  the  surface 
to  an  amorphous  condition.  ^  Under  such  conditions  it  was 
found  that  the  laser  beam  melts  the  crystal  to  a  depth  some¬ 
what  greater  than  that  of  the  implanted  profile.  This  molten 
layer  recrystallizes  from  the  undamaged  substrate  generating 
an  epitaxially  grown  defect-free  single  crystalline  silicon 

layer  containing  the  implanted  impurity  atoms  in  normal 
( 5 ) 

lattice  sites.  Thus  both  defect  elimination  and  electrical 

activation  of  the  dopant  atoms  are  simultaneously  carried 
out.  During  the  time  the  implanted  region  is  molten,  the 
impurity  atoms  diffuse  rapidly  in  the  liquid,  resulting  in  a 
significant  redistribution  of  the  implanted  dopant  atoms. 


It  has  been  suggested  that  the  time  available  for  annealing 
with  Q-switched  ruby  lasers  is  too  short  for  an  unmolten 
damaged  layer  to  be  annealed  with  a  single  pulse  and  that 
melting  followed  by  liquid  phase  epitaxy  is  essential  for 
effective  annealing. ^ 

Continuous  Ar  laser  annealing  was  applied  by  A.  Gat 

et  al.  ^  to  silicon  surfaces  implanted  with  fluences  of 

14  -2 

As  heavy  enough, 5x10  cm  at  100  keV,  to  convert  the  surface 
to  an  amorphous  condition.  Melting  did  not  occur  and  thus  also 
no  redistribution  of  the  implanted  As.  Fewer  defects  remained 
than  were  found  for  a  single  step  thermal  anneal  (30  min  at 
1000°C)  but  complete  elimination  of  defect  nuclei  did  not 
occur . 

14 

In  this  study  we  chose  B  implant  fluences  of  10  , 

14  15  -2  14 

5x10  and  10  cm  and  P  implant  fluences  of  10  and 

14  -2 

5x10  cm  into  both  {100}  and  {111}  silicon  surfaces.  One 

14  -2 

implantation,  that  of  5x10  Pcm  into  a  {111}  surface, 
resulted  in  amorphism.  For  the  rest  of  the  wafers,  the 
surfaces  retained  their  single  crystallinity. 

For  these  lower  fluence  implants  it  was  found  that 
laser  energies  which  were  described  in  the  literature  as 
resulting  in  melting  of  the  surface  and  completely  activating 
the  implanted  impurity  atoms  here  failed  to  heat  the  silicon 
surfaces  sufficiently  to  permit  more  than  a  fraction  of  the 


implanted  dopant  atoms  to  migrate  to  substitutional  positions. 
Nevertheless  the  temperature  distribution  within  the  implanted 
layer  was  such  that  complete  solution  of  implant  defect  nuclei 
occurred.  By  combining  a  laser  anneal  with  a  short  thermal 
anneal  it  was  possible  to  achieve  both  defect  elimination  as 
well  as  complete  electrical  activation.  Such  a  process  would 
avoid  impurity  redistribution,  permitting  the  dopant  distri¬ 
bution  to  remain  as  implanted. 

II.  EXPERIMENTAL  INVESTIGATION 

Twenty  wafers  were  subjected  to  implantation  with  boron 
and  phosphorus  ions.  The  orientations,  conductivity  types  and 
resistivities  of  the  wafers, together  with  the  fluences  and 
implant  energies  are  given  in  Table  I.  The  implant  energies 
were  chosen  to  give  similar  distributions  for  the  B  and  P 
ions.  The  wafers  were  commercially  polished  and  had  resis¬ 
tivities  in  the  range  of  1-4  Dcm.  After  implantation  the 
wafers  were  examined  for  the  presence  of  amorphism.  This  was 
determined  by  the  presence  of  a  characteristic  "milky" 

appearance  of  the  implanted  area  as  compared  to  the  adjacent 

(2) 

areas  protected  from  implantation  by  the  wafer  holder. 

Wafers  216A  and  216B  exhibited  such  amorphism. 

After  implantation  the  series  A  wafers  were  subjected 
to  an  identical  laser  annealing  process  and  the  series  B 
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wafers  were  halved,  each  half  being  subjected  to  a  separate 
laser  annealing  process.  Figure  1  describes  the  experimental 
procedure . 

The  laser  annealing  was  carried  out  with  a  Q-switched 

ruby  laser  operated  in  the  multimode  condition.  The  pulse 

duration  was  15  nsecs  and  the  beam  size  was  approximately 

1  cm  in  diameter.  The  beam  was  projected  through  a  rectangular 

aperture  4.4  mm  by  7.9  mm.  For  the  A  wafers,  large  areas, 

approximately  one  fourth  of  the  surface  of  each  2  inch  diameter 

wafer,  were  pulsed  with  a  single  laser  energy.  This  was 

accomplished  by  using  pulses  of  the  same  energy, as  the  wafer 

was  repeatedly  translated  a  distance  equal  to  the  width  or 

height  of  the  projected  beam.  Figure  2  illustrates  the  laser 

pulse  distribution.  The  laser  energies  utilized  were  1.0,  1.5, 
-2 

2.0  and  2.5  J  cm  .  After  laser  annealing , all  surfaces  were 
examined  by  high  power  optical  microscopy.  The  rippled 
appearance  characteristic  of  surface  melting  ^  was  not 
found.  In  the  center  of  each  wafer  there  remained  a  small 
reference  area  which  was  not  subjected  to  laser  annealing. 

The  laser  annealed  areas  were  sufficiently  large  to  permit 
a  sheet  resistance  value  to  be  obtained  directly  without 
introducing  correction  factors. 

Tables  II  and  III  list  the  results  of  the  sheet  resis¬ 
tance  measurements  following  the  laser  annealing.  A  minimum 
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of  6  readings  was  taken  for  each  annealed  area  and  the  average 
taken.  For  the  two  lower  energies  there  was  considerable 
scatter  while  for  the  higher  two  energies  all  readings  were 
within  +  5%  of  the  mean.  This  would  appear  possibly  to  be 
due  to  inhomogeneities  within  the  laser  beam.  For  the  higher 
energies,  lateral  thermal  conduction  could  smooth  out  the 
annealing  effect. 

After  measuring  the  sheet  resistances  of  the  laser 
annealed  areas  all  wafers  were  subjected  to  a  thermal  anneal 
of  10  minutes  at  1000°C  in  an  N2  atmosphere.  The  results 
of  sheet  resistance  measurements  repeated  on  the  same  laser 
annealed  areas  are  found  in  Table  IV. 

As  seen  in  Fig.  1,  the  B  series  wafers  were  halved, 

one  half  designated  as  B1  and  the  other  as  B2.  The  Bl  halves 

were  subjected  to  laser  annealing  at  9  discrete  areas.  Each 

area  was  4.4  mm  x  7.9  mm,  the  size  of  the  aperture  through 

which  the  beam  was  projected.  The  laser  energies  used  were 

_  o 

.5,  1.0,  1.5,  2.0  and  2.5  J  cm  with  a  single  pulse  as  well 
as  .5  (2X) ,  .5  (4X) ,  1.0  (2X) ,  and  1.0  J  cm2  (3X) .  Each 
discrete  laser  annealed  area  was  surrounded  by  a  non  annealed 
area  and  after  oxidation  and  Sirtl  etching,  the  laser  pulse 
boundary  (LPB)  was  examined  by  optical  microscopy  and  compared 
to  the  ion  implant  boundary  (IIB)  between  the  area  of  the  wafer 
covered  by  the  wafer  holder  and  the  area  of  the  wafer  subjected 
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to  implantation  without  laser  annealing.  The  results  are  listed 
in  Table  V. 

Wafers  116B  and  216B  received  identical  implantations 
14  -2 

of  5x10  Pcm  ,  but  only  216B  with  a  {111  }  orientation  exhibited 
amorphism.  The  wafer  halves  116B2  and  216B2  were  subjected  to  a 
pre-anneal  treatment  of  1  h  at  550°C  to  permit  slow  trans¬ 
formation  of  the  amorphous  layer  to  single  crystallinity.  This 
was  confirmed  by  the  disappearance  of  the  "milky"  surface  on 
216B2.  All  B2  halves  were  then  annealed  for  10  min  at  1000°C. 
Sheet  resistivity  measurements  were  carried  out  and  are  given 
in  Table  IV  under  the  heading  of  zero  laser  beam  energy. 

The  B2  halves  were  then  subjected  to  laser  annealing  at 

5  discrete  areas,  using  single  laser  pulsed  of  .5,  1.0,  1.5, 

_  o 

2.0  and  2.5  J  cm  .  Each  discrete  laser  annealed  area, 

4.4  mm  x  7.9  mm,  was  surrounded  by  a  non  annealed  area.  After 
steam  oxidation  and  Sirtl  etching,  the  surfaces  were  examined 
by  optical  microscopy.  The  results  are  listed  in  Table  VI. 

III.  EXPERIMENTAL  RESULTS 

Examination  of  the  results  listed  in  Tables  II  and  III 

shows  that  in  every  case  there  is  a  significant  increase  in 

electrical  activity  of  the  implanted  ions  in  going  from  1.0 
—  2  —  2 

to  1.5  J  cm  and  from  1.5  to  2.0  J  cm  .  In  going  from  2.0 

__  O 

to  2.5  J  cm  ,  the  change  in  Rp  appears  to  be  much  smaller. 
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At  first  it  might  appear  that  this  is  a  result  of  an  approach 
to  complete  activation.  Following  the  laser  anneal  with  the 
moderate  thermal  anneal  of  10  minutes  at  1000°C  results  in 
a  very  significant  decrease  in  sheet  resistance  for  all  speci 
mens  for  all  laser  anneals  as  is  shown  by  comparing  Table  IV 
with  Tables  II  and  III.  This  is  illustrated  in  Fig.  3,  in 
which  the  sheet  resistances  for  laser  annealing  alone  and 
for  combined  laser  and  thermal  annealing  are  compared.  We 
recognize  for  the  fluences  studied  that  laser  annealing  by 
itself  is  incapable  of  bringing  about  full  electrical  acti¬ 
vation  of  the  implanted  ions. 

A  measure  of  the  effectiveness  of  laser  pulses  in 

eliminating  implant  damage  is  the  comparison  of  the  implanted 

and  laser  annealed  surfaces  with  that  of  an  unimplanted  and 

unpulsed  surface  after  oxidation  and  Sirtl  etching.  Where 

these  are  identical  we  can  reasonably  assume  that  the  laser 

annealing  eliminated  the  defect  nuclei  introduced  by  the  ion 

implantation.  As  seen  in  Table  V,  this  occurred  for  most  of 

-2 

the  B1  wafer  halves  for  laser  pulses  of  1.5  J  cm  and  for  all 

-2 

of  them  for  2.0  and  2.5  J  cm  .  For  a  laser  pulse  of  1.0  J  cm 
there  was  a  transitional  structure  between  no  annealing  and 
complete  annealing.  This  structure  was  also  found  at  the 
LPB  of  the  higher  energy  laser  annealed  areas. 


Figure  4  illustrates  the  laser  pulse  boundary  (LPB) 

_  2 

for  specimen  113B1  for  1.5  J  cm  .  The  non-pulsed  surface 

is  characterized  by  a  light  concentration  of  oxidation 

stacking  faults  (OSF) ,  the  pulsed  area  is  relatively  clear, 

and  the  transitional  region  is  characterized  by  a  higher 

concentration  of  OSF  corresponding  to  that  found  for  the 

_o 

lower  laser  energy  of  1.0  J  cm  .  Figures  5  and  6  illustrate 

_  o 

the  LPB  for  specimen  117B1  for  1.5  J  cm  .  In  Fig.  5  we 

see  the  laser  annealed  area  free  of  OSF  on  the  left  and  the 

transition  zone  containing  a  heavy  concentration  of  OSF  with 

a  light  background  on  the  right.  In  Fig.  6  the  transition 

zone  is  on  the  left  and  the  unannealed  area  containing  a 

heavy  background  together  with  a  heavy  concentration  of  OSF 

—  2 

is  found  on  the  right.  For  a  laser  energy  of  1.0  J  cm  ,  the 

annealed  area  is  very  similar  to  the  transition  area  for  the 
-2 

1.5  J  cm  anneal.  No  results  are  listed  for  the  lowest  laser 

-2 

energy,. 5  J  cm  ,  since  no  differences  in  ternary  defect 
distribution  from  the  non  pulsed  areas  was  observed. 

In  this  investigation  we  also  studied  the  effect  of 

-2 

multiple  laser  pulses  of  .5  and  1.0  J  cm  .  From  Table  V, 

_2 

we  note  that  a  laser  pulse  of  1.0  J  cm  results  in  a  signifi¬ 
cant  reduction  of  the  defect  concentration  for  various  im- 

-2 

plantations  while  1.5  J  cm  was  effective  in  eliminating 
all  defects  in  most  cases.  It  appeared  reasonable  that 
multiple  pulses  of  1.0  J  cm  might  be  equivalent  to  a  single 
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pulse  of  higher  energy.  It  was  found  however,  that  multiple 

-2 

pulses  of  .5  and  1.0  J  cm  did  not  reduce  the  defect  con¬ 
centration  to  any  significant  extent  below  the  level  found  for 
a  single  pulse  of  the  same  energy. 

Following  the  thermal  anneal  of  the  A  series  wafers, 
these  coo  were  subjected  to  steam  oxidation  and  Sirtl 
etching.  Microscopic  examination  showed  results  similar  to 
those  obtained  with  the  Bl  series,  confirming  our  prior  ex¬ 
perience  that  an  anneal  of  1000°C  for  10  minutes  has  little 
or  no  effect  in  dissolving  ion  implant  defect  nuclei. 

The  annealing  procedure,  10  min  at  1000°C  in  dry  N^, 
utilized  for  the  B2  wafer  halves  prior  to  laser  pulsing 
converts  the  ion  implant  damage  into  secondary  defects.  ^ 
These  consist  of  small  dislocation  loops,  probably  resulting 
from  the  agglomeration  of  the  primary  defects,  the  vacancies 
and  interstitial  atoms  generated  in  the  implantation  process. 
The  secondary  defects  appear  to  be  much  more  stable  than  the 
primary  defects  as  shown  by  their  greater  resistance  to 
solution  by  laser  annealing.  These  results  are  listed  in 
Table  VI.  For  three  of  the  lower  fluence  specimens,  113B2, 
213B2,  and  214B2,  the  expanded  implant  defect  concentration 
was  not  sufficiently  greater  than  the  background  to  permit 

us  to  detect  clearly  a  LPB.  Comparing  Table  V  with  Table  VI 

-2 

we  note  that, where  in  the  former  a  laser  energy  of  2.0  J  cm 
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was  sufficient  to  dissolve  the  defect  nuclei  in  every  single 

specimen,  in  the  latter  we  find  that  the  higher  energy  of 
-2 

2.5  J  cm  results  only  in  partial  reduction  of  the  defect 
concentration  with  one  notable  exception,  specimen  216B2. 

Wafers  216A  and  B  were  the  only  ones  of  the  20  studied 
which  exhibited  amorphism  after  ion  implantation.  Specimens 
116B2  and  216B2  received  identical  two  step  thermal  anneals 
prior  to  laser  pulsing.  The  first  step.  60  min  at  550°C,  was 
chosen  to  permit  slow  regrowth  of  the  amorphous  layer  into 
single  crystal  silicon.  Examination  of  specimen  216B2  after 
this  step  did  indeed  show  that  the  implanted  surface  had  lost 
its  "milky"  appearance  and  could  not  be  distinguished  from 
the  adjacent  non- implanted  surfaces.  Sheet  resistivity 
measurements  gave  values  of  392  n/r,  for  116B2  and  264  fi/n 
for  216B2,  agreeing  with  work  of  Miyao  et  al v  that  only 
partial  activation  is  obtained  for  550°C  anneals. 

The  second  step,  10  min  at  1000°C,  was  a  standard 
thermal  anneal  used  to  produce  electrical  activation  of  the 
implanted  impurities.  The  two  step  process  has  been  shown  to 
reduce  significantly  the  implant  damage  for  heavier  implant 

(9) 

fluences 

We  found  that  the  use  of  a  two  step  thermal  anneal  for 

this  amorphous  surface  permits  us  to  dissolve  the  ion  implant 

-2 

defect  nuclei  with  a  laser  pulse  energy  of  2.0  J  cm  . 


This  last  result  suggests  that  the  technique  of  combined 
laser  and  thermal  annealing  which  was  developed  for  low 
fluence  implants  as  described  here  might  also  be  effectively 
applied  to  heavy  implants  where  the  surface  is  converted  to  an 
amorphous  condition  and  therefore  prone  to  melting  and  con¬ 
comitant  impurity  redistribution  when  subjected  to  pulsed 
laser  annealing.  By  introducing  a  two  step  thermal  anneal, 
the  amorphous  surface  is  solid  phase  epitaxially  regrown 
without  impurity  redistribution.  Subsequent  laser  annealing 
results  in  the  dissolution  of  the  implant  defects  without 
melting  or  the  resultant  impurity  redistribution. 


IV.  DISCUSSION 

The  results  of  this  study  suggest  that  the  energy  release 
in  the  interaction  of  a  pulsed  ruby  laser  with  a  single 
crystalline  surface  containing  discrete  implant  defects  occurs 
primarily  at  the  defects,  creating  localized  hot  spots  which 
result  in  the  dissolution  of  the  defect  and  in  the  heating  and 
electrical  activation  of  the  area  adjacent  to  the  defect.  This 
phenomenon  was  clearly  illustrated  in  a  separate  unpublished 
study  of  the  interaction  of  Q  switch  ruby  laser  pulses, 
identical  to  the  ones  used  in  this  study, with  silicon  which 
had  been  subjected  to  neutron  bombardment.  In  that  study 
we  found  clear  signs  of  the  occurrence  of  localized  melting 


at  the  silicon  surfaces  in  the  vicinity  of  defects.  Fig.  7. 
We  did  not  detect  such  melt  spots  for  ion  implanted  and 
laser  annealed  silicon. 

The  interaction  of  a  laser  pulse  with  a  semiconducting 
crystal  is  well  understood.  For  each  photon  absorbed  we 
generate  a  hole-electron  pair.  These  have  a  kinetic  energy 
as  well  as  a  potential  energy  equal  to  the  band  gap  energy. 
The  kinetic  energy  is  transferred  almost  immediately  to  the 
crystal  lattice  as  heat,  while  the  potential  energy  is  re¬ 
leased  at  such  a  time  and  place  as  recombination  of  the 
hole-electron  pair  occurs. 

We  recognize  that  the  intensity  of  a  laser  pulse,  I, 
in  penetrating  a  distance  x  from  the  surface  of  a  homogene¬ 
ous  crystal  is  given  by  the  relation 

_  _  -  a  x 

I  =  I  e 
o 

where 

Iq  is  the  initial  intensity  and 

a  is  the  absorptivity. 

The  concentration  of  excess  hole-electron  pairs 
generated  is  proportional  to  the  rate  of  absorption 

M  =  O-  =  '  ii  =  aI  (x) 

The  surface  volume  contains  not  only  the  implanted 
impurities  but  also  a  high  concentration  of  crystal  defects 
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which  can  act  as  recombination  centers  for  the  hole-electron 
pairs  generated  by  the  laser  pulse.  This  permits  the  gener¬ 
ation  of  high  temperatures  at  defect  sites  by  the  diffusion 
of  excess  electrons  and  holes  from  adjacent  areas.  A  require¬ 
ment  for  this  model  to  be  applicable  is  that  the  rate  of 
diffusion  of  holes  and  electrons  to  the  defect  recombination 
centers  exceed  by  a  significant  amount  the  rate  of  thermal 
diffusion  from  the  defect. 

Another  model  views  the  defect  site  as  important  in  the 
generation  of  hole-electron  pairs. 

The  absorptivity,  a,  is  a  function  of  the  crystalline 
perfection  of  the  silicon.  The  value  for  amorphous  silicon 
exceeds  that  for  crystalline  silicon  by  an  order  of  magnitude. 
This  suggests  that  the  absorption  of  laser  energy  at  a  defect 
site  is  much  higher  than  in  adjacent  undamaged  silicon.  Ac¬ 
cording  to  this  view,  we  would  generate  high  concentrations 
of  holes  and  electrons  around  individual  defect  sites.  We 
recognize  that  the  recombination  rate  is  proportional  to 
the  hole-electron  concentration  product,  i.e.,  R  =  r  (^nj^p^j  . 
An  increase  in  the  electron  and  hole  concentrations  of  an 
order  of  magnitude  results  in  an  increase  in  the  recombination 
rate  of  two  orders  of  magnitude.  The  rate  of  diffusion  of 
excess  holes  and  electrons  away  from  the  defect  site  where 
they  were  generated  will  be  proportional  to  their  concentration 
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gradient  and  hence  to  their  concentration.  This  is  signifi¬ 
cantly  less  than  the  rate  of  recombination,  insuring  that 
recombination  and  hence  thermal  release  occurs  at  the  defect 
sites  where  the  hole-electron  pairs  were  generated. 

On  the  basis  of  available  evidence  we  are  unable  to 
decide  whether  either  or  both  of  these  suggested  models  are 
valid. 

In  our  results  we  found  that  repeated  lower  energy 

-2 

pulses,  i.e.,  .5  and  1.0  J  cm  resulted  in  no  greater  solution 
of  defect  nuclei  than  occurs  with  a  single  pulse.  It  would 
appear  that  the  first  pulse  dissolves  some  of  the  primary 
defect  nuclei  while  converting  others  due  to  released  thermal 
energy  into  secondary  defects.  As  shown  in  our  second  set 
of  experiments,  secondary  defects  are  much  more  resistant 
to  dissolution  by  laser  pulsing,  so  that  repeated  pulsing 
contributes  little  to  defect  removal. 

V.  CONCLUSIONS 

1.  Pulsed  laser  annealing  of  lower  fluence  ion  implants 

does  not  result  in  full  electrical  activation. 

-2 

2.  Laser  anneal  single  pulse  energies  of  1.5  J  cm  are 
sufficient  to  dissolve  implant  defect  nuclei  in 
lower  fluence  ion  implants  without  causing  melting. 

3.  Following  a  pulsed  laser  anneal  with  thermal 
annealing  permits  us  to  obtain  both  dissolution 
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of  implant  defect  nuclei  as  well  as  complete  electrical 
activation  without  redistribution  of  low  fluence  ion 
implanted  impurities. 

4.  A  two  step  thermal  anneal  followed  by  a  pulsed  laser 
anneal  permits  us  to  obtain  both  dissolution  of  implant 
defect  nuclei,  as  well  as  complete  electrical  activation 
without  impurity  redistribution, for  implant  fluences 
sufficiently  high  to  cause  amorphism  at  the  surface. 

This  work  was  supported  in  part  by  the  Defense  Advanced  Research 
Projects  Agency. 
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TABLE  I 

EXPERIMENTAL  CONDITIONS 


WAFER  NUMBERS  TYPE  ORIENTATION 
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SPECIES 


IMPLANT 
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Boron  implantations  at  50  keV 
Phosphorus  implantations  at  80  keV 


TABLE  II 


SHEET  RESISTANCES  OF  BORON  IMPLANTED 
WAFERS  FOLLOWING  LASER  ANNEAL 


SHEET  RESISTANCE  fl/D 
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NUMBER 
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TABLE  III 

SHEET  RESISTANCES  OF  PHOSPHORUS  IMPLANTED 
WAFERS  FOLLOWING  LASER  ANNEAL 


SHEET  RESISTANCE.  (l/C 
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SHEET  RESISTANCE  IN  fl/Q  OF  LASER  ANNEALED  AREAS 
FOLLOWING  THERMAL  ANNEAL  OF  10  MINUTES 
AT  1000°C  IN  N2  AMBIENT 
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Figure  1.  Flow  Diagram  of  Experiments 


Figure  2.  Distribution  of  Laser  Pulses.  Area  I  (1.0 
Area  II  (1.5  J  cm' 2),  Area  III  (2.0  J  cm'2 
Area  IV  (2.5  J  an* 2). 


LASER  BEAM  ENERGY  (J  cm  ) 


Figure  3 


Canparison  of  Laser  Annealing  and  Laser  Plus  Thermal 
Annealing  of  Sheet  Resistances  of  5  x  10^  Bern' 2 
Implanted  Wafers 


Figure  4. 


Laser  Pulse  Boundary  for  Specimen  113B1. 
lx  10  B  cm  %  1.5  J  on'-  Laser  Energy 
Sirtl  Etch  10  Secs,  Mag.  100X. 
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Figure  5. 


Laser  Pulse  Boundary  for  Specimen  Uni¬ 
junction  of  Laser  .Annealed  and  Transition 
Area.  1  x  1015  B  cm'2,  1.5  J  cm'-  Laser 
Energy,  Sirtl  Etch  10  Secs,  Mag.  250X. 
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Figure  6. 


Laser  Pulse  Boundary  for  Specimen  117B1. 
Junction,  Unannealed,  and  Transition  Area. 

c-X  ?n  ’  1,5  J  cm’“  Laser  Energy, 
Sirtl  Etch  10  Secs,  Mag.  250X. 
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Figure  7.  Scanning  Electron  Micrograph  of  Neutron 
Irradiated  Silicon  Surface  Exposed  to  a 
foiby  Laser  Pulse  of  2.0  J  cm"-,  Mag.  300X. 
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